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Abstract 

Aims. We aim to determine the fundamental parameters of a sample of B stars with apparent visual magnitudes below 8 in the field-of-view of 
, the CoRoT space mission, from high-resolution spectroscopy. 

Methods. We developed an automatic procedure for the spectroscopic analysis of B-type stars with winds, based on an extensive grid of 
I FASTWIND model atmospheres. We use the equivalent widths and/or the line profile shapes of continuum normalized hydrogen, helium and 
D silicon line profiles to determine the fundamental properties of these stars in an automated way. 

' Results. After thorough tests, both on synthetic datasets and on very high-quality, high-resolution spectra of B stars for which we already had 
, accurate values of their physical properties from alternative analyses, we applied our method to 66 B-type stars contained in the ground-based 
I archive of the CoRoT space mission. We discuss the statistical properties of the sample and compare them with those predicted by evolutionary 
' models of B stars. 

, Conclusions. Our spectroscopic results provide a valuable starting point for any future seismic modelling of the stars, should they be observed 
■ by CoRoT. 

Key words. Stars: atmospheres - Stars: early-type - Stars: fundamental parameters - Methods: data analysis - Techniques: spectroscopic - 
' Line: profiles 



1. Introduction 

The detailed spectroscopic analysis of B-type stars has for 
a long time been restricted to a limited number of targets. 
Reasons for this are the a priori need for a realistic atmosphere 
model, the lack of large samples with high-quality spectra, and 
the long-winded process of line profile fitting, as the multitude 
of photospheric and wind parameters requires a large parameter 
space to be explored. 

The advent of high-resolution, high signal-to-noise spec- 
troscopy in the nineties led to a renewed interest of the sci- 
entific community in spectroscopic research, and in partic- 
ular in the relatively poorly understood massive stars. The 
establishment of continuously better instrumentation and the 
improvement in quality of the obtained spectroscopic data 
triggered a series of studies, which led to a rapid increase 
in our knowledge of massive stars. In this respect, it is not 
surprising to note that this is exactly the period where sev- 
eral groups started to upgrade their atmosphere prediction 



code for such stars, see, e.g., CMFGEN - Hillier & Millet 
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PHOENIX - 



Hauschildt & Baron' ('1999'), WM-Basic 



Pauldrach et all (I2001h. POWR - .Grafener et al.. (.2002) an d 



FASTWIND - ISantolava-Rev et aLlTiggTl) . IPuIs et al]l2005b . 
Initially, major attention was devoted to the establishment of 
a realistic atmosphere model (improvement of atomic data, in- 
clusion of line blanketing and clumping), rather than analyzing 
large samples of stars. 

Simultaneously with improvements in the atmosphere pre- 
dictions, also the number of available high-quality data in- 
creased rapidly, mainly thanks to the advent of multi-object 
spectroscopy. At the time of writing, the l argest survey is th e 
VLT-FLAMES Survey of Massive Stars dEvans et al.l l2005h . 
containing over 600 Galactic, SMC and LMC B-type spectra' 
(in 7 different clusters), gathered over more than 100 hours 
of VLT time. The survey not only allowed to derive the stel- 
lar parameters and rotational velocities for hundreds of stars 
dOufton et al.ll2006t iHunter et al]|2008bh . but also to study the 
evolution of surface N abundances and the effective tempera- 



additionally, roughly 90 0-stars have been observed. 
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ture scales in the Galaxy and Ma gellanic Clouds (ITrundle et al 
2007tlHunter et al.ll2007il2008ah . 



In preparation of the CoRoT space mission, and almost 
contemporary with the FLAMES setup, another large database 
was constructed: GAUDI (G round-based Asteroseismology 



Uniform Database Interface, ISolano et alJ l2005h . It gathers 



ground-based observations of more than 1500 objects, includ- 
ing high-resolution spectra of about 250 massive B-type stars, 
with the goal to determine their fundamental parameters as in- 
put for seismic modeling (see Section[3]). 

The availability of such large samples of B-type stars brings 
within reach different types of studies, e.g., they may lead to 
a significant improvement in the fundamental parameter cali- 
bration for this temperature range and to a confr ontation with 



and ev aluation of stellar evolution models (e.g., iHunter et al 



2008bh . The drawback of this huge flood of data, however, is. 



as mentioned before, the large parameter space to be explored, 
which can be quite time-consuming, if no adequate method is 
available. It requires a method which is able to derive the com- 
plete set of parameters of stars with a wide variety of physical 
properties in an objective way. 

To deal with the large GAUDI dataset, we investigated the 
possibility of automated spectral line fitting and we opted for 
a grid-based fitting method; AnalyseBstar. In SectionlH we 
justify our choice for a grid-based method, present its design 
and discuss several tests which were applied to check the per- 
formance of the routine. A more detailed description of our 
methodology can be found in the (online) appendix. Section[3] 
illustrates the first application of AnalyseBstar to the sample of 
CoRoT candidate targets in the GAUDI database and Section|4] 
deals with the physical interpretation and some statistical prop- 
erties of the resulting parameters. Section|5] summarizes the 
main results obtained in this paper. 

2. Automated fitting using a grid-based methiod 

Spectral line fitting is a clear example of an optimization prob- 
lem. To find the optimal fit to a given observed stellar spec- 
trum among a set of theoretically predicted spectra emerging 
from stellar atmosphere models, requires scanning the param- 
eter space spanned by the free parameters of the stellar atmo- 
sphere model. Due to the extent of this parameter space, and 
with the goal of analyzing large samples in mind, it is clear 
that performing the parameter scan through fit-by-eye is not the 
best option. The 'subjective' eye should be replaced by an 'in- 
telligent' algorithm, in such way that the procedure of finding 
the optimum fit becomes automatic, objective, fast and repro- 
ducible, even though human intervention can never be excluded 
completely. 



2. 1. General description of our grid-based mettiod 

In contrast to the case of the O stars (e.g.. lMokiem et al .l2005h . 
treating the entire spectral range B with a single approach re- 
quires additional diagnostic lines besides H and He. In this 
region. Si (in its different ionization stages), rather than He, 
becomes the most appropriate temperature indicator With this 
in mind, we have chosen to develop an automatic procedure 



which is based on an extensive and refined grid of FASTWIND 
models. This offers a good compromise between effort, time 
and precision if an appropriate grid has been set up. The grid 
should be comprehensive, as dense as possible and representa- 
tive for the kind of objects one wants to analyze. Similar to a 
fit-by-eye method, a grid-based algorithm will follow an iter- 
ative scheme, but in a reproducible way using a goodness-of- 
fit parameter. Starting from a first guess for the fundamental 
parameters, based on spectral type and/or published informa- 
tion, improved solutions are derived by comparing line profiles 
resulting from well-chosen existing (i.e. pre-calculated) grid 
models to the observed line profiles. The algorithm terminates 
once the fit quality cannot be improved anymore by modifying 
the model parameters. 

We are well aware of the limitations inherent to this 
method. As soon as a new version/update of the atmosphere 
or line synthesis code is released (e.g., due to improved atomic 
data), the grid needs to be updated as well. The advantage, on 
the other hand, is that the line-profile fitting method itself will 
remain and the job is done with the computation of a new grid. 
This is of course a huge work (e.g., seven months were ne eded 



to compute and check the grid in iLefever et al.l 12007 al) . but 



thanks to the fast performance of the FASTWIND code, this 
should not really be an insurmountable problem. 

A grid-based method fully relies on a static grid and no ad- 
ditional models are computed to derive the most likely param- 
eters. Consequently, the quality of the final best fit and the pre- 
cision of the final physical parameters are fully determined by 
the density of the grid. This underlines, again, the necessity of 
a grid which is as dense as possible, allowing for interpolation. 
On the other hand, the grid-method has a plus-point: it is fast. 
Whereas the analysis with, e.g., a genetic algorithm approach 
(IMokiem et al.ll2005l e.g.) needs parallel processors and sev- 
eral days of CPU to treat one target star, our grid-method will, 
on average, require less than half an hour on one computer, be- 
cause no additional model computations are required. 

2.2. Tine code AnalyseBstar 

Our grid-based code, called AnalyseBstar, was developed for 
the spectroscopic analysis of B stars with winds. It is written 
in the Interactive Data Language (IDL), which allows for inter- 
active manipulation and visualization of data. It fully relies on 
the extensive grid of NLTE m odel atmosphe r es and t he emerg- 
ing line profiles presented in Lefever et al. ( 2007a[ see also 
Appendix lB . SI for a brief overview of the considered param- 
eter set) and is developed to treat large samples of stars in a 
homogeneous way. We use continuum normalized H, He and 
Si lines to derive the photospheric properties of the star and 
the characteristics of its wind. Such an automated method is 



- 180 CPU months were needed to compute the full grid of al- 
most 265 000 models. To reduce the effective computation time, cal- 
culations were done on a dedicated Linux cluster of 5 dual-core, 
dual-processor computers (3800 MHz processors, sharing 4 Gb RAM 
memory and 8 Gb swap memory), amounting to 20 dedicated CPU 
processors, in addition to 40 more regularly used institute CPUs (8 
of 3800MHz and 32 of 3400MHz). The grid filled 60% of a terabyte 
disk, connected to a Solaris 10 host pc. 
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not only homogeneous, but also objective and robust. A full 
description of the design of AnalyseBstar, including the prepa- 
ration of the input, the iteration cycle for the determination of 
each physical parameter and the inherent assumptions, is given 
in Appendix B-*. 

The intrinsic nature of our procedure allows us to derive es- 
timates for the 'real' parameters of the star, which in most cases 
are located in between two model grid points. In what follows, 
we refer to them as 'interpolated values'. The corresponding 
grid values are the parameters corresponding to the grid model 
lying closest to these interpolated values. 

Furthermore, we will denote the surface gravity as derived 
from fitting the Balmer line wings as \ogg, while the gravity 
corrected for centrifugal terms (required, e.g., to estimate con- 
sistent masses) will be denoted by \og gc- The (approximate) 
correction itself is obtain ed by adding the term (vsin / )^/Rt to 
the uncorrected gravity (iRepolust et al.ll2005l and references 
therein). 

2.3. Testing the method 

2.3.1 . Convergence tests for synthetic FASTWIND 
spectra 

Before applying AnalyseBstar to observed stellar spectra, we 
tested whether the method was able to recover the parameters 
of synthetic FASTWIND input spectra. To this end, we created 
several synthetic datasets in various regions of parameter space, 
with properties representative for a 'typical' GAUDI spectrum. 
For more details on the setup of this test dataset and the re- 
sults of the convergence tests, we refer to Appendix|A] In all 
cases, the input parameters were well recovered and no conver- 
gence problems were encountered. Minor deviations from the 
input parameters were as expected and within the error bars. 
Also the derived v sin / -values, which are difficult to disentan- 
gle from potential macroturbulent velocities, agreed very well 
with the inserted values, irrespecti ve of Vmacro- This sh ows that 
the Fourier Transf orm method of iGrav d 1 97311 1975k . whose 
implementation bv lSimon-Dfaz & Herrero (l2007h we used, in- 
deed allows to separate both effects (but see below). All to- 
gether, this gave us enough confidence to believe that our pro- 
cedure will also be able to recover the physical parameters from 
real spectral data. 

2.3.2. FASTWIND vs. Tlusty: profile comparison 
and analysis of Tlusty spectra 

As an additional test, we compared line profiles based on 
FASTWIND models with negligible mass loss with their cor- 
responding counterpart in the grid computed with the NLTE 
model atmosphere code Tlusty (iHubenv & Lanz l2000l) . which 
allows for a fully consistent NLTE metal line blanketing as 
well. This resulted in an overall good agreement, except for 
low gravity stars, where differences arise. It concerns differ- 
ences in the forbidden He I components (which are quire strong 
in FASTWIND, but almost absent in Tlusty), the Balmer line 



wings (which are a bit stronger in FASTWIND than in Tlusty, 
yielding slightly lower logg values) and the Sill/lll EW ra- 
tio for Teff > 18 kK (which is larger for FASTWIND than 
for Tlusty - due to both stronger Si II and weaker Si III lines -, 
yielding higher reft values). The reason that FASTWIND pre- 
dicts stronger Balmer line wings is due to differences in the 
broadening functions and an underestimated photospheric line 
pressure, which can affect the photospheric structure of low- 
gravity objects. Differences are about or below 0.1 dex, and 
become negligible for dwarf stars. Except for the problem with 
the forbidden components, the profiles of the cooler models 
(Tfiff < 18 kK) compare very well. 

To check the uncertainty in the synthetic profiles, we ap- 
plied AnalyseBstar to profiles of six prototypical Tlusty models 
at 3 different temperature points (15 000, 20 000, and 25 000 K) 
and at a high and a low gravity. We added artificial noise to the 
data and artificially broadened the profiles with a fixed pro- 
jected rotational velocity. The results of the comparison can be 
found in Table[T] The input and output fundamental parameters 
agree very well within the resulting errors, with only one ex- 
ception, being the model with the high Tf^g and higher logg, 
where the temperature is off by 2 000 K. 

Around 20kK, we experience similar difficulties to retrieve 
a unique temperature as in the cool B star domain, i.e. depend- 
ing on the gravity, we may have only one stage of Si available, 
in this case Si III. Indeed, whereas the Si II lines are clearly vis- 
ible in the higher gravity test case (logg = 3.00), they are no 
longer detectable in the lower gravity case (logg = 2.25). At 
this Teff, there are (unfortunately) no obvious Si IV lines yet, 
which complicates the analysis. Even though, AnalyseBstar is 
still able to retrieve the correct effective temperature, albeit 
with somewhat larger errors, using the same alternative method 
as for late B type stars when there is only Si II available (see 
'method 2' in Section lBTJl l. 



Table 1. Result of applying AnalyseBstar to synthetic Tlusty 
profiles, which are artificially adapted to obtain a typical SNR 
of 150 and a projected rotational broadening of vsin/ = 50 
kms Teft.in and log gin are, respectively, the input effec- 
tive temperature and gravity of the Tlusty spectra, reft^out and 
log g out represent the best matching (Teg^, log g)-combination 
in the FASTWIND grid. 



T'eff.in 




7'cfF.out ± A 


logg.out ± A 


(kK) 


(cgs) 


(kK) 


(cgs) 


15 


1.75 


15.0 ±0.5 


1.7 ±0.05 


15 


3.00 


15.5 ±0.5 


3.0 ±0.10 


20 


2.25 


20.0 ± 1.0 


2.2 ±0.10 


20 


3.00 


21.0 ± 1.0 


3.2 ± 0.20 


25 


2.75 


25.0 ± 1.0 


2.7 ± 0.05 


25 


3.00 


27.0 ± 1.0 


3.1 ±0.10 



^ Manual available upon request from the authors. 
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15 20 25 30 2.5 3.0 3.5 4.0 4.5 

(kK) - DETAIL/SURFACE log g (DETAIL/SUREACE) 



Figure 1. Comparis on of the effective tempera tures (left) and surface gravities (right) derived with FASTWIND (AnalyseBstar) 
with those found bv lMorel etal] ( l2006l |2008l) and iBriquet & Morell OOP?) using DETAIL/SURFACE, for the jS Cephei stars 
(black), SPBs (grey) and some well-studied hot B stars (open symbols). Supergiants are indicated as diamonds, giants as triangles 
and dwarfs as circles. Th e shaded ar ea around the one-to- one relation ( straight line) denotes the uncertainty on the derived 
fundamental parameters of iMorel et al., C2006. .2008.) and .Briquet & Morell (l2007i) (i.e., 1 000 K in Teff and 0.15 in \ogg). 



2.3.3. Comparison with fit-by-eye results for 
well-studied pulsators 

Having thoroughly tested that AnalyseBstar indeed converges 
towards the optimum solution and that it is able to recover the 
input parameters of synthetic spectra, we performed an ad- 
ditional test on real spectra. We tested our method on a se- 
lected sample of high-quality, high-resolution spectra of pulsat- 
ing B stars (fi Cephei and Slowly Pulsating B stars (SPB)). The 
(mean) spectra of these stars have a very high SNR, attained 
thro ugh the addition of a large number of individual exposures 
(see. iMorel et al.ll2007l) . Given that we excluded binaries from 
our GAUDI sample (see Section lTTI ). we also ignored 6 Oph 
and j8 Cru for the test sample. Moreover, V2052 Oph is chemi- 
cally peculiar and was thus excluded, as well as V836 Cen for 
which only one spectrum is available. Thus we limited our test 
sample to those single stars in Morel et al. (2007) for which 
numerous high-quality spectra were available. Because of the 
high quality of these spectral time series, they are ideally suited 
for testing AnalyseBstar. 

The p Ce phei stars an d two of the S PB st ars were analyzed 
in det ail by iMorel etaP (|2006[ l2008h and iBriquet & Morell 
(l2007h . respectively. These authors used the latest version 
of the NLTE line formation codes DETAIL and SURFACE 



dGiddingsl 119811: iButler & Giddingsl Il985h . in combination 



with plane-parallel, fully l ine-blanketed LTE Kurucz atmo- 
spheric models (ATLAS9, iKuruczl 119931) . to determine (by 
eye) the atmospheric parameters and element abundances 
of low-luminosity class objects with negligible winds. The 
outcome of the comparison with their results are summa- 
rized in Tablesl2al to |2c] and in Figs.[T] and |2] Although 
the effective temperatures derived from FASTWIND tend to 
be, on the average, slightly below the ones derived from 
DETAIL/SURFACE, it is clear that, within the error bars, 
FASTWIND and DETAIL/SURFACE give consistent results 
for Teff and \ogg (Fig.[T]) as well as for log n(Si)/n(H) (Fig.|2]i, 



-4.2 




-5.0 -4.8 -4.6 -4.4 -4.2 -4.0 
log n(Si)/n(H) (DETAIL/SURFACE) 



3.8 



Figure 2. Comparison between the Si abundances 
derived from AnalyseBstar and those derived from 
DETAIL/SURFACE for the set of photometric targets, used to 
evaluate the performance of AnalyseBstar The black dotted 
lines represent the solar Si abundance (log n(Si)/n(H) = -4.49), 
while the grey dotted lines represent a depletion and enhance- 
ment of Si by 0.3 dex (i.e., -4.79 and -4.19, respectively). 
Symbols have the same meaning as in Fig.[T] Individual error 
bars are indicated. 



On th e other hand, the v sin / values derived by IMorel et al 



(120081) are typically slightly above those derived by us. This is 
readily understood as their values implicitly include the macro- 
turbulent velocity, which has not been taken into account as a 
separate broadening component. In those cases where we de- 
rive a vanishing V r naCTn, t he v sin / values are in perfect agree- 



ment. 



Aerts et al.l ( 120091) recently suggested that macroturbu- 



lence might be explained in terms of collective pulsational 
velocity broadening due to the superposition of a multitude 
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Table 2a. Fundamental parameters of well-studied p Cephei stars (spectral types from SIMBAD), as derived from AnalyseBstar: 
the effective temperature {T^g), the surface gravity {\ogg, not corrected for centrifugal acceleration due to stellar rotation), the 
helium abundance (n(He)/n(H), with solar value 0.10), the Si abundance (log n(Si)/n(H), with solar value -4.49), the microtur- 
bulence (^), the projected rotational velocity (vsin/ ) and the macroturbulent velocity (Vmacro)- We list the values of the clos- 
est grid model, except for the Si abundance, where the in terpolated value is given. In italics, we display the results from the 
DETAIL/SURFACE analysis by iMorel et alJ (|2006[ 120081) . For 12 Lac, we additionally list the parameters which we retrieve 
when forcing Vmacio to be zero while keeping the value for vsin ; (indicated between brackets, see text for details). All stars have 
thin or negligible winds. 

Typical errors for reff an d \ogg are, respectively, 1 000 K and 0.15 for the comparison data sets of iMorel et al. I (l2006ll2008l) and 
Briquet & Morell (120071) . which is somewhat more conservative than the errors adopted in this study. We adopt typical 1-cr errors 
of 0.10 dex for log g, 1 000 K for > 20 000 K and Teff < 15 000 K, and 500 K for 15 000 K < T^g < 20 000 K. The error 
for the derived Si-abundance, log n(Si)/n(H), has been estimated as 0.15 and 0.20 dex for objects above and below 15 000 K, 
respectively (see also Section lB3] ) 



HD 


alternative 


Spectral 


T^ff 




n(He)/n(H) 


log n(Si)/n(H) 




V sin i 




number 


name 


Type 


(K) 


(cgs) 








(km s-') 


(kms 


-') 


(kms-') 


46328 


f ' CMa 


B0.5-B1IV 


27000 


3.80 





10 


-4.69 


6 


9± 


2 


11 








27500 


3.75 






-4.87 ±0.21 


6±2 


10 ± 


2 




50707 


15 CMa 


BlIII 


24000 


3.40 





10 


-4.79 


12 


34 ± 


4 


38 








26000 


3.60 






-4.69 ± 0.30 


7±3 


45 ± 


3 




205 021 


ySCep 


BlIV 


25000 


3.80 





10 


-4.70 


6 


26 ± 


3 


24 








26000 


3.70 






-4.89 ± 0.23 


6±3 


29 ± 


2 




44743 


ySCMa 


B1.5III 


24000 


3.50 





10 


-4.76 


15 


19 ± 


4 


20 








24000 


3.50 






-4.83 ± 0.23 


14 ±3 


23 ± 


2 




214993 


12 Lac 


B1.5IV 


23000 


3.60 





10 


-4.41 


6 


44 ± 


6 


37 
















[-4.83] 


[12] 






[0] 








24500 


3.65 






-4.89 ± 0.27 


10 ±4 


42 ± 


4 




16582 


(JCeti 


B1.5-B2IV 


23000 


3.90 





10 


-4.80 


6 


14 ± 


2 











23000 


3.80 






-4.72 ± 0.29 




14 ± 


1 




886 


7 Peg 


B1.5-B2IV 


23000 


3.80 





10 


-4.84 


<3 


10 ± 


1 











22500 


3.75 






-4.81 ± 0.29 




10 ± 


1 




29248 


y Eri 


B1.5-B2IV 


23000 


3.70 





10 


-4.73 


10 


21 ± 


3 


39 








23500 


3.75 






-4.79 + 0.26 


10 ±4 


36 ± 


3 





Table 2b. Same as for Table | 2al bu t now for some well-studied SPBs. In italics, we show the results of the DETAIL/SURFACE 
analysis by iBriquet & IVIorell (l2007l) . 



HD 


alternative 


Spectral 




log,? 


n(He)/n(H) 


log n(Si)/n(H) 




V sin i 


^macro 


number 


name 


Type 


(K) 


(cgs) 








(kms"') 


(km s-') 


(kms-') 


3360 


^ Cas 


B2IV 


22000 


3.80 





10 


-4.76 


<3 


18 + 2 


13 








22000 


3.70 






-4.72 ±0.30 


1 ± 1 


19 ±1 




85953 


V335 Vel 


B2IV 


20000 


3.80 





10 


-4.84 


6 


30 ± 1 


20 








21000 


3.80 






-4.75 ±0.30 


1±1 


29 ±2 




3379 


53 Psc 


B2.5IV 


20000 


4.30 





10 


-4.73 


<3 


48 ±8 


43 


74195 


Vel 


B3IV 


16500 


3.70 





10 


-4.79 


<3 


18 ± 2 


18 


160762 


I Her 


B3IV 


19500 


4.10 





10 


-4.86 


3 


8± 2 





25558 


40Tau 


B3V 


17500 


4.00 





10 


-4.72 


6 


28 ± 2 


31 


181558 


HR 7339 


B5III 


15000 


4.00 





10 


-4.79 


6 


17 ± 2 





26326 


HR 1288 


B5IV 


15500 


3.60 





10 


-4.49 


<3 


17 ± 1 


17 


206540 


HR 8292 


B5IV 


13500 


3.80 





10 


-4.79 


3 


15 ± 2 





24587 


HR 1213 


B5V 


14500 


4.00 





10 


-4.79 


<3 


25 ±4 


21 


28114 


HR 1397 


B6IV 


14000 


3.50 





10 


-4.79 


<3 


21 ±4 


17 


138764 


HR 5780 


B6IV 


14500 


3.90 





10 


-4.51 


3 


21 ± 2 





215573 


HR 8663 


B6IV 


13500 


3.80 





10 


-4.49 


<3 


8± 1 





39844 


HR 2064 


B6V 


14500 


3.70 





10 


-4.39 


<3 


16 ± 1 





191295 


V1473 Aql 


B7III 


13000 


3.70 





10 


-4.79 


3 


16 ± 1 


15 


21071 


V576 Per 


B7V 


13500 


3.70 





10 


-4.79 


3 


22 ± 1 





37151 


VI 179 Ori 


B8V 


12500 


3.80 





10 


-4.79 


<3 


20 ±2 
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Table 2c. Same as for Table|2al but now for some wel l -studi ed hot B stars. In italics, we display the independen t results from 
the DE TAIL/SURFACE analysis by iMorel et al] ilOOd I2OO8I) . excep t for 6 Car and t Sco which can be found in lHubrig etaT 



(I2OO8I) . For r Sco, we additionally compare with the values found by iMokiem et alj (12005 ) usin g a genetic algorithm approach 
(superscript M). The results for two objects overlapping with the study of Przybilla et alj ( 2008 ) are also indicated (superscript 
P). 



HD 


alternative 


Spectral 




log.? 


n(He)/n(H) 


log n(Si)/n(H) 






V sin / 


^macro 


number 


name 


Type 


(K) 


(cgs) 






(km s 




(km s-') 


(km s-') 


93030 


6/ Car 


BOVp 


31000 


4.20 


0.15 


-4.38 




10 


108 ±3 











31000 


4.20 




-4.57 ± 0.23 


12 : 


t4 


113 ± 8 




149 438 


T Sco 


B0.2V 


32000 


4.20 


0.10 


-4.55 




6 


10 ±2 











31500 


4.05 




-4.76 ±0.14 


2: 


t2 


8±2 










31900" 


4.15" 


0.12^ 




10. 


gM 












32000'' 


4.3(f 


0.10'' 


-4.50'' 




5'' 


4'' 


4'' 


36591 


HR 1861 


BIV 


26000 


3.90 


0.10 


-4.73 




6 


13 ± 1 











27000 


4.00 




-4.75 ±0.29 


3z 


t2 


16 ±2 










27000'' 


4.12'' 


O.W 


-4.52'' 




3'' 


W 




52089 


6CMa 


B1.5-B2II/III 


22000 


3.20 


0.10 


-4.69 




15 


32 ±2 


20 








23000 


3.30 




-4.77 ±0.24 


16 z 


b4 


28 ±2 




35468 


7 Ori 


B2n-in 


22000 


3.60 


0.10 


-4.79 




10 


46 ±8 


37 








22000 


3.50 




-5.00 ±0.19 


13 z 


t5 


51 ±4 




51309 


(CMa 


B2.5Ib-II 


17000 


2.60 


0.10 


-4.78 




15 


27 ±4 


39 








17500 


2.75 




-4.82 ± 0.31 


15 z 


t5 


32 ±3 





of gravity modes with low amplitudes (see also Appendix 
IB ■ 41 ). In this respect, we expect stars for which Vmacro differs 
from zero to be pulsators. 

For few overlapping objects, we could als o compare our 



results with those from independent studies by iMokiem et al 



(120051) . based on FASTWIND predictions for H and He hnes 



and us ing a genetic algorithm approach, and by iPrzvbilla et al 
(I2OO8I) . using ATLAS9/DETAIL/SURFACE, allowing for a 



comparison of H, He, and Si. The consistency with both stud- 
ies is quite good, except for the Si abundance of HD 36591 



as dete rmined by Przybilla et al. For this object, iHubrig et al 



(l2008h have obtained a value close to ours, using the same 
ATLAS9/DETAIL/SURFACE code. Ad hoc, we cannot judge 
the origin of this discrepancy. 

For two targets, HD 85953 and t Sco, somewhat larger 
differences in the derived microturbulence are found. For 
HD 85953, this difference is compensated by an opposite 
difference in the derived Si abundance. For r Sco, on the 
other hand, the derived microturbulence (and also the Si abun- 
dance) is in perfect agreement with th e values provided by 
Hubrig et all (120081) and Przybilla et al.l (l2008i) within the de- 



fined error bars, whereas the value for the microturbulence de- 



rived by IMokiem et al.l (120051) differs substantially. As they did 
not derive the Si abundance, we cannot judge on the difference 
in Si abundance this would have caused. 

An even more detailed test can be done for stars with seis- 
micaUy determined values of the fundamental parameters. 12 
Lac is such a star It is known as a non-radial pulsator with 



at lea st 1 1 independent oscillation frequencies ([Handler et al 



20061) . The fact that 12 Lac is a rich pulsator clearly shows up 



in the line profiles, which are strongly asymmetric. Especially 
the Si III triplets (Si III 4552-4567-4574 and Si III 4813-4819- 
4829) are skew. Fig. [3] shows the observed spectrum (black), 
which is a combination of 31 exposures. Using the Fourier 



Transform method, we find the projected rotational velocity to 
be 44 kms"', which is slightly above th e 36 kms~^ derive d 
from modelling the fine profile variations (iDesmet et al ]l2009h . 
Including macroturbulence in the spectral line fitting, results in 
the following parameters: T^g = 23 000 K ± 1 000 K, log g = 3.6 
±0.1,^ = 6 + 3 km s ' , solar He and Si abundances (solid grey 
lines in Fig.[3]i. Actually, the derived surface gravity is in per- 
fect agreement with the s eismically determine d value of log g 
(3.64-3.70) determined bv lDesmet eTal] (|2009|) . 



Excep t for the Si abunda nce, our results are also in agree- 
ment with lMorel et al.l (120061 see Tablel2ali. who derived Teff = 
24500K± 1 000 K,logg = 3.65 + 0.15,^= 10 + 4 kms-' and 
depleted Si abundances as seem to be typical for B dwarfs in 
the solar neighborhood. To account for the pulsational broad- 
ening, we needed to include a macroturbulence of 37 kms 
Forcing the macroturbulence to be zero while keeping the de- 
rived value of vsin / leaves the main physical parameters (Teff, 
\ogg and He abundance) unaltered, while the microturbulent 
velocity and the Si abundance change to compensate for the 
change in profile shape: the microturbulence becomes larger 
(12 km s"', broader lines) and the Si abundance becomes lower 
(depleted, log n(Si)/n(H) a -4.83). As expected, t his de pleted 
Si abundance is in agreement with iMorel et al.l (120061) . who 
did not include Vmacro in their analysis. When evaluating the 
fit quality, we find, besides the expected mismatch in the Si 
line wings, also a discrepancy i n the line cores of He (dashed 
grey lines in Fig.|3]l. Following lAerts et al.l (120091) . this might 
be improved by accounting for the pulsational broadening in- 
troduced by the collective set of detected oscillations in 12 Lac, 
as this would result in a shape intermediate between a rotation 
and a Gaussian profile. However, the only way one can account 
for this, is by using an appropriate time series of spectra. The 
above example illustrates that we need to account for pulsa- 
tional broadening to correctly assess the Si abundance and mi- 
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cro turbulence. In this sense, it is advisable to work with ' Vmacio' 
as a substitute for the time-dependent broadening in the case of 
'snapshot' spectra. 

3. Application to thie GAUDI B star sample 

3.1. C0R0T& GAUDI 

The French-led European space mission CoRoT (Convection, 
Rotation, and plane tarv Transits; see "the CoRoT Book", 
Fridlund et aDl2006l) was launched successfully on December 
27th, 2006. The observational setup of the CoRoT seismol- 
ogy programme (observations of a small number of bright 
stars over a very long period) made target selection a cru- 
cial issue. Due to the acute shortage of available information 
on the potential CoRoT targets, the need for additional data 
was high. Therefore, an ambitious ground-based observing pro- 
gram for more than 1500 objects was set up, under the leader- 
ship of C. Catala at Meudon, to obtain Stromgren photome- 
try {uvbyp) as well as high-resolution spectroscopy (FEROS, 
ELODIE, SARG, CORALIE, GIRAFFE, coude spectrograph 
at the 2m telescope in Tautenburg, CATANIA). These data 
were collected in an extensive catalogue, maintained at LAEFF 
(Laboratorio de Astrofisica Espacial y Fisica Fundamental) and 
baptized GAUDI: Ground-based Aste roseismology Uniform 
Database Interface (I Solano et al ] l2005h . 

The Be s t ars wi thin t he GAUDI database were studied by 
Fremat et al. ( 20061) and Neiner et al. (2005) and are omitted 
here, since the FASTWIND code is not developed to treat stars 
with a circumstellar disc. Also double-lined spectroscopic bi- 
naries (SB2) were omitted from our sample, as their combined 
spectra make an accurate fundamental parameter estimation 
impossible as long as their flux ratios are not known. Among 
the GAUDI sample, we discovered a few candidate spectro- 
scopic binaries, which were not known to be SB2s. It con- 
cerns HD 173003 (B5), HD 181474 (B5), HD 42959 (B8), 
HD 50982 (B8, variable star), HD 51 150 (B8, clearly asymmet- 
ric profiles), HD 181761 (B8), HD 45953 (B9, variable star) 
and HD 46165 (B9). 

For the few available SARG data, only the default normal- 
ized spectra were inserted into the database and we failed to 
get hands on the raw data, which would be needed for a care- 
ful rerectification. Indeed, the poor continuum rectification is 
especially clear from H/?, but also shows up in other line pro- 
files. Moreover, the spectral coverage of the SARG spectra is 
too small to deduce any useful information. Therefore, we ex- 
cluded all SARG spectra from our sample. The same accounts 
for the CATANIA spectra, for which the spectral quality was 
too poor for a detailed spectral analysis, so also these spec- 
tra were omitted. We thus restrict this paper to FEROS and 
ELODIE spectra, having a resolution of 48 000 and 50 000, re- 
spectively. 

The standard FITS data of the FEROS and ELODIE spectra 
in the GAUDI database contain information about both the nor- 
malized spectrum, resulting from the pipeline reduction, and 
the unnormalized spectrum. Since the quality of normalization 
turned out to be insufficient for our detailed analyses, one of 
us (TM) redid the normalization for all spectral ranges around 



the diagnostic lines in a uniform way using IRAF"* to ensure a 
homogeneous treatment of the sample. 

Among the GAUDI sample, we find, as expected, a lot 
of fast rotators. Their spectra usually contain insufficient line 
information due to their high projected rotational velocities. 
Often, only the Balmer lines and the strongest He I lines (i.e. 
He I 4026, 4471 and 4922) can be detected, while all Si 
Unes are lost. The stars that were discarded for this reason 
are HD 182519 (B5, >300 kms '), HD 50751 (B8, >200 
kms-i), HD 56006 (B8, >150kms-i), HD 45515 (B8 V, >190 
kms '), HD 182786 (B8, >160 kms '), HD 50252 (B9 V, 
>140kms-'), HD 169225 (B9, >160 kms '), HD 171931 (B9, 
>250 kms-'), HD 174836 (B9, >140 kms"'), 176258 (B9 V, 
> 160 km s '), HD 179124 (B9 V, >270 km s '), and HD 45760 
(B9.5 V, >190kms '). For some of the GAUDI stars, we were 
able to fit the spectrum despite the high projected rotational ve- 
locity, but their resulting parameters cannot be very reliable. 
This concerns HD 51507 (B3 V, 148 kms"'), HD 45418 (B5, 
237 kms-i), HD 46487 (B5 Vn, 265 kms"'), HD 178744 (B5 
Vn, 224 kms-'), HD 43461 (B6 V, 210 kms"'), HD 44720 
(B8, 160 km s '), HD 49643 (B8 Illn, 296 km s '), HD 173370 
(B9V, 282 km s ', double-peaked Ha profile), HD 179124 (B9 
V, 278 km s '), and HD 181690 (B9 V, 189 kms '). Their pa- 
rameters can be found in Table|3] 

For other targets, we have only spectra of insufficient qual- 
ity available which cannot be used for spectral line fitting pur- 
poses. It concerns HD 48691 (B0.5 IV), HD 168797 (B3 Ve), 
HD 178129 (B3 la), HD 57608 (B8 III, possible instrumen- 
tal problem), HD 44654 (B9), HD 45257 (B9), and HD 53204 
(B9). 

For a few stars, we have both an ELODIE and a FEROS 
spectrum available, e.g., for HD 174069. From inspection of 
the spectra, it was immediately clear that there are differences 
in the line profiles, in particular the wings of the Balmer lines 
are much less pronounced in the ELODIE spectra (see Fig.|4|l. 
Also other stars, for which we have both a FEROS and an 
ELODIE spectrum available, show the same discrepancy, so 
there seems to be a systematic effect. We fitted both spec- 
tra and came up with a different set of parameters. The dif- 
ference in \ogg is large (> 0.5 dex) and certainly worrisome. 
To investigate this problem, we looked up the spectra in the 
FEROS and ELODIE archives and realized that the merging of 
the ELODIE spectra by the pipeline is far less accurate than 
for the FEROS spectra. This is due to an inaccurate correc- 
tion for the blaze function in the case of the ELODIE pipeline. 
Due to this, we also redid the merging of the orders for the 
ELODIE spectra, before the normalization. This led to a much 
better agreement with the available FEROS spectra. We thus 
advise future users against working with the merged and nor- 
malized spectra from GAUDI, but rather to go back to the orig- 
inal spectra in the FEROS and ELODIE archives, and not only 
redo the normalization but even the merging in the case of the 
ELODIE spectra. Unfortunately, not all stars had a spectrum 



* IRAF (Image Reduction and Analysis Facility) is distributed 
by the National Optical Astronomy Observatories, operated by the 
Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation, USA. 
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Silll4574 




4549.2 4551.4 4553.6 4555.7 
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4572.9 4574.1 4575.4 4576.6 



HALPHA 




HBETA 




HEI4922 




Silll481 3 



6515 6546 



4821 4848 4875 4902 



4918 4926 




11.5 4812.7 4813.9 4815.2 



Silll4819 




17.8 4819.1 4820.4 4821.7 



SIIII4829 




4826.6 4828.1 4829.6 4831.1 



HDELTA 




4064 4089 4115 4140 



1 .00 

0, 



Sill41 28 



4126.67 4127.51 4128.56 4129.50 



Sill4130 




4129.46 4130.44 4131.41 4132.38 



SilV41 1 6 




4114.3 41 15.4 4116.5 4117.7 



HGAMMA 




4307 4329 4352 4375 



HEI6678 




71.4 6675.7 6680.0 6684.3 




5737.2 5738.8 5740.5 5742.1 



0.94 b . . 

5039 5040 5042 5043 





4379.6 4385.1 4390.7 4396.2 



0, 
0, 
0, 
0, 

4014.4 4021.8 4029.3 4036.7 



Figure 3. Illustration of the effects of time-independent broadening on the derived parameters for 12 Lac (B2 III). The observed 
line profiles (solid black lines) are fitted using AnalyseBstar Inclusion of time-independent broadening (by means of macrotur- 
bulence) yields the solid grey fit, with optimum values of v sin ; = 44 km s"' and Vmacio = 37 km s"' . Forcing the macroturbulence 
to be zero leaves all parameters but the Si abundance and the microturbulence unaltered. The resulting fit is represented by the 
grey dashed lines. 



available in the ELODIE archives though, e.g., for HD 47887 
(B2 III), HD 52559 (B2 IV- V), HD 48977 (B2.5 V), HD 50228 
(B5), HD 57291 (B5), HD 51892 (B7 III), HD 52206 (B8) and 
HD 53202 (B9), we could not find the original spectrum, so we 
left out these stars from our analysis. For HD 43317 (B3 IV), 
the archival spectrum was not usable. 

Finally, for 17 of the remaining stars, we could not find a 
satisfactory fit, and we decided to leave them out of the sample. 
One possible explanation might be that they are single-lined 
binaries (see also iMassey et al., 2009; who encountered simi- 
lar problems when analyzing a large sample of LMC/SMC O- 
stars). 

All together, this, unfortunately, reduces our sample a lot, 
and we are only left with a bit more than a third of the targets 
initially in the database: 66 out of 187 objects. Moreover, most 
of the objects are late B-type stars. Fig.|5] shows a histogram 
with the breakdown of the sample following spectral subtype. 
The ratio marked above each bin indicates the number of ana- 
lyzed targets compared to the total amount that was available 
before the selection procedure. 



3.2. Analysis results 
3.2.1. Example fits 

To illustrate the obtained quality of the final fits, we show in 
Figs.|6]to[8]some examples of the resulting spectral fits for three 
very different stars included in GAUDI: a 'cool' giant, a hot 
supergiant and a middle type dwarf. 
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4840 4850 4860 4870 4880 
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Figure 4. HD 174069; Significant (and apparently system- 
atic) discrepancies are observed between the line profiles of 
the FEROS (solid black line) and the ELODIE (dashed black 
line) spectrum from the GAUDI database. This leads to consid- 
erable discrepancies in the derived stellar parameters. In grey, 
we show the normalized ELODIE spectrum obtained from the 
ELODIE archives, which shows a much better agreement with 
the available FEROS spectra. 



10 



Lefever et al.: Fundamental parameters of B stars in the FoV of CoRoT 



Table 3: Stellar parameters of the analyzed GAUDI B type stars. The first two columns give the HD number of the star and the spectral type (SpT) taken from SIMBAD. Then we 
list the effective temperature (T^tr) of the model which was found by AnalyseBstar to match the observations best, the surface gravity corrected for centrifugal acceleration (log^^, 
with coiTesponding uncertainty Alog^}. the He abundance (n(He}/n(H), where 0. 10 is solar), the interpolated Si abundance when available, otherwise the closest grid Si abundance 
(log n(Si)/n(H), where -4.49 is solar, -4.19 is enha nced and -4 . 79 is dep leted) and the microturbulenc e (^); vsin/ is the projected rotational velocity as obtained by applying the 
implementation of the Fourier Transform method of I Gray! 119731 n'97^ bv lSim6i> Diaz & Herrerol |2Q07|) . In some cases, we need an additional broadening to explain the line profiles. 
Therefore, we list the required macroturbulent velocity (Vmacro) and indicate the stars for which we suspect they might be pulsators from the fact that they have a significant VmLiao- 
We attribute a flag to each star, to represent the obtained fit quality or problems we encountered during the fitting procedure: 1 - well-fitted, 2 - probably a wrong spectral type (about 
1 000 K or more away from the calibrations, proposed spectral type indicated in the 'remarks'-column). 3 - no or hardly any Si lines available, mostly due to fast rotation, 4 - not very 
reliable fit (for other reasons than fast rotation, see 'remarks'), 5 - He weak/Si strong, 6 - fitting of Ha, Yifi, and Hy line profiles complicated due to very broad wings which abruptly 
change into very narrow line cores (log g may be less accurate). In the last column, a few additional comments are given. The formal uncertainty on the effective temperature is 1 ,000 K 
for ToiT > 20, 000 K and T^n < 1 5 000 K, and 500 K for T^n < 1 5, 000 K < T^^x < 20 000 K. The formal uncertainty on the microturbulence and on the He and Si abundance is half of 
the grid step in most cases: 2 km s"' and 0.025 dex, respectively. For objects below 15 000 K we adopt somewhat largers enors (3-iT deviation con'esponding to 2 grid steps): 3 km s"\ 
0.033 and 0.2 dex, respectively. 



HD number SpT T^{f log gc n(He)/n(H) log n(Si)/n(H) ^ v sin / Vmucro suspected flag remarks 

(K) (cgs) (km s~') (km s~') (km s~^ ) pulsator 



48434 


BOIII 


28000 


3.11 ±0.10 


O.IO 


-4.44 


15 


62 ±3 


31 


yes 


4 


ELODIE, poor quality spectrum. He lines too strong 
Balmer lines not well fitted due to bumps in blue wings 


52382 


Bllb 


23000 


2.71 ± 0.10 


0.15 


-4.67 


20 


56 ±5 


53 


yes 


4 


P Cygni profile Ha not reproduced, 
H/? and Hy not well fitted 


170580 


B2V 


20000 


4.10 ± 0.10 


0.10 


-4.81 


6 


11 ± I 





no 


I 


N(He)/N(H) lower than solar 


44700 


B3V 


17500 


3.80 ±0.10 


O.IO 


-4.93 


6 


8± 1 





no 


I 


N(He)/N(H) lower than solar 


181074 


B3 


20000 


3.64 ±0.10 


O.IO 


-4.35 


12 


133 ±7 





no 


2 


probably B2, fast rotator 


45418 


B5 


16000 


4.0 ± 0.20 


O.IO 


-4.49 


3 


237 ± 28 





no 


3 




46487 


B5Vn 


14500 


3.69 ±0.10 


0.15 


-4.49 


3 


265 ±5 





no 


3 




48215 


B5V 


14500 


3.82 ±0.10 


O.IO 


-4.49 


3 


98 ±3 


48 


yes 


1 


only Si 11 


54596 


B5 


20000 


3.31 ±0.05 


0.10 


-4.48 


6 


69 ±6 


53 


yes 


2 


only Si 111, probably B2 


58973 


B5 


15000 


3.43 ± 0.20 


0.10 


-4.49 


3 


94 ± 1 


50 


yes 


1 




177880 


B5V 


14500 


3.81 ±0.10 


0.10 


-4.49 


3 


49 ± 1 


23 


yes 


1 




178744 


B5Vn 


14500 


3.71 ±0.20 


0.10 


-4.49 


12 


224 ±5 





no 


3 




43461 


B6V 


13500 


3.35 ±0.10 


0.10 


-4.79 


3 


210 ±5 





no 


3 




48807 


B7Iab 


12500 


2.00 ±0.10 


0.10 


-4.28 


6 


24 ± 1 


27 


yes 


1 




51360 


B71I1 


13500 


3.12 ±0.10 


0.10 


-3.93 


3 


73 ± 1 


55 


yes 


1 


very high Si abundance 


42677 


B8 


11000 


3.55 ±0.20 


0.10 


-4.79 


6 


121 ± 14 


104 


yes 


2 


probably B9 


44720 


B8 


14000 


3.96 ±0.10 


0.10 


-4.19 


3 


160 ±5 





no 


2 


probably B5 or B6 


45153 


B8 


1 1000 


3.66 ±0.10 


0.10 


-4.49 


3 


142 ±4 





no 


2 


probably B9 


45284 


B8 


15000 


4. 10 ±0.10 


0.10 


-4.19 


3 


52 ±5 


42 


yes 


2 


SPB with magnetic field, skew Si 11 lines, probably B5 


45397 


B8 


11500 


3.59 ±0.10 


0.20 


-4.79 


6 


156 ± 12 





no 


3 




45515 


B8V 


11000 


4.18 ±0.10 


0.15 


-4.49 


3 


190 ±5 





no 


2, 3 


probably B9 


46616 


B8 


17500 


4.20 ± 0.50 


0.10 


-4.24 


3 


8± 1 


7 


yes 


5 


He-weak star 


47964 


B8II1 


11500 


3.01 ± 0.20 


0.10 


-4.79 


3 


49 ± 1 


44 


yes 


1 




48497 


B8 


14000 


3.70 ±0.10 


O.IO 


-4.58 


3 


14 ± 1 


13 


yes 


4,6 


N(He)/N(H) lower than solar 


49481 


B8 


11000 


2.70 ±0.10 


0.20 


-4.79 


3 


9± 1 


13 


yes 


2,6 


probably B9 


49643 


BSlIln 


14500 


3.88 ±0.10 


0.10 


-4.19 


3 


296 ±2 


13 


no 


3 




49886 


B8 


10000 


3.30 ±0.20 


O.IO 


-4.79 


6 


8± 1 


16 


yes 


4,6 




49935 


B8 


12500 


3.33 ± 0.20 


0.10 


-4.79 


6 


92 ±9 


72 


yes 


1 


log g may be too low? 


50251 


B8V 


11500 


3.00 ± 0.20 


0.10 


-4.79 


3 


9± 1 


17 


yes 


1 


emission line around Si 111 4813 and Si IV 4212 


50513 


B8 


11500 


4.04 ±0.10 


0.10 


-4.79 


6 


112± 5 


80 


yes 


1 




55793 


B8 


11000 


3.05 ± 0.05 


0.20 


-4.79 


3 


101 ± 1 


56 


yes 


2 


probably B9, log^ may be too low? 


56446 


B81I1 


11500 


3.21 ±0.20 


0.20 


-4.19 


3 


229 ± 32 


13 


no 


3 




170795 


B8 


15500 


4.01 ±0.10 


0.10 


-4.49 


3 


79 ±7 


37 


yes 


2 


probably B5 


171247 


B81Ilsp 


10000 


2.83 ±0.10 


0.10 


-4.19 


6 


66 ±2 


20 


yes 


4 


peculiar line behaviour 


173673 


B8 


10000 


2.90 ±0.10 


0.20 


-4.49 


3 


25 ±2 


26 


yes 


4,6 


too low log 


179761 


B811-111 


12500 


3.30 ± 0.05 


0.10 


-4.49 


3 


17 ± 1 


12 


yes 


I 
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B8 


17000 
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0.10 
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6 


167 ±3 





no 


3 


wrong spectral type? 
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30 
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I 
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0.10 
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12 


no 
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no 
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10 


no 
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B9 
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56 


yes 


I 
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3 
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no 


4 
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B9 
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3 
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23 


yes 


1 


log g too low? 


47278 


B9 
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0.10 
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3 


31 ±4 


25 


yes 


2 


rather AO 
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B9 
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0.10 


-4.79 


3 


47 ± 14 


63 


yes 


1 
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3 
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23 
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3 
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no 


2 


rather AO 
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3 
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no 


3 
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54 
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B9 
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6 


53 ± 1 


71 


yes 


1 


Balmer line cores slightly refilled, log g too high? 
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B9 
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0.10 


-4.49 


3 


46 ± 1 


107 


yes 


4 


fast rotator, difficult to fit, too high Vmacro 


56613 


B9 
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O.IO 
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6 


91 ±7 


36 


yes 


2 


rather, B7 or B8 
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B9 
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3.62 ±0.10 


0.10 
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3 
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no 


1 




173370 


B9V 


11500 


3.46 ±0.10 


0.10 


-4.79 


15 


282 ±8 





no 


4 


very fast rotator, no Si lines available,double-peaked Ho- profile, 
H/5 refilled on both sides of the line core 


173693 


B9 
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3 
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22 
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rather AO 
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3 
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3 
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3 
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yes 


2 


rather AO 
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3 
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no 


2 


rather AO 
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3 
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no 


2 


rather B 8 


179124 


B9V 


12000 


3.45 ±0.10 
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no 


3 
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32 
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1 
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13000 


3.42 ±0.10 


0.10 
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6 
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no 


3 
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B9V 


11000 


3.10 ±0.10 


O.IO 
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3 


23 ± 1 


12 


yes 


I 
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Spectral type (taken from SIMBAD) 

Figure 5. Histogram of the analyzed stars within the GAUDI 
B star sample. The ratio marked above each bin indicates the 
number of analyzed targets compared to the total amount that 
was available before the selection procedure, e.g., 1/5 means 
that there were 5 objects with this spectral type in the database, 
but that, for certain reasons, we were only able to analyze 1 of 
them. In total, 66 out of 187 objects could be analyzed. 



The corresponding physical parameters are listed in 
Table[3] which summarizes the derived stellar parameters of 
the analysed GAUDI B type stars. All other spectral Une 
fits, and their corresponding parameters, can be found at 
|http://www.ster.kuleuven.be/~karohen/AnalyseBstar/Bstars/ 

3.2.2. Comments on individual stars 

In what follows, we discuss a few of the sample stars individu- 
ally, because of certain particularities, 

HD 45284. The derived effective t emperature is in com plete 
agreement with what was found by iHubrig et al. I (I2OO6I) . This 
star is an SPB exhibiting a magnetic field. The combination 
of the effecti ve temperature and gravity resulting from Geneva 
photometry dHubrig et al. I I2OO6I) positions the star just outside 
the main sequence towards the high-gravity side. Now, with 
our newly derived log g, it falls within the expected instability 
domain of the SPB stars. 



HD 4661 6 (B8) is a clear example of a He-weak star: all He I- 
lines have extremely low equivalent widths compared to 'nor- 
mal' B stars. Even though the hydrogen and Si lines fit nicely, 
the He lines are predicted too strong in our grid models^. 

HD 481 06 (B8) is oxygen poor based on the non-detection of 
the O I 7771-7775 triplet and is either He weak or Si strong. 
We were unable to fit both the He and Si lines at the same time. 

For HD 46340 (B8), only He 1 4471 and Si 11 5041-5056 are 
visible. This was only sufficient to make a very rough estimate 
of the fundamental parameters: T^ff ^ 10000 K, \ogg ^ 4.1, 
negligible wind, solar He abundance. 



^ the lowermost considered He-abundance in our grid is 
n(He)/n(H)=0.1, see SectionlB^Sl 



On the hot side of the B-type domain, the analysis of the 
early type stars (e.g., HD 172488 (B0.5 V), HD 52918 (Bl V) 
and HD 173198 (Bl V)) show similar difficulties. The iterative 
procedure tends to yield too low temperatures due to the lower 
weight of the Si IV and He II lines as we have only one weak 
line. As only Si III is reliable and He I is not much affected by 
the effective temperature in this range, we are not able to derive 
trustworthy results. 



HD 1 73370 (B9 V) shows a double peaked Ho- profile. This 
morphology is typical for a very fast rotator (in this case, 
vsin/ » 280 kms '), seen almost equator-on. The disk-fike 
feature also strongly affects the H/3 and Hy profiles. 

We also found a few chemically peculiar stars among the 
B stars in the GAUDI database. They are HD 44948 (B8 Vp), 
HD 45583 (B8), HD 46616 (B8) and HD 44907 (B9). They 
have a very rich spectrum, with a forest of sharp spectral lines, 
complicating the continuum determination and the fitting pro- 
cess, as almost all lines are blended. The He lines of HD 45583 
are extremely weak and those of HD 46616 even completely 
vanished. Following Landstreet (2007), HD 45583 is a period- 
ically variable Ap star The line spectrum is completely dis- 
torted, hence too complicated to fit. The Sill lines are very 
strong, which indeed points towards a B8 or even B9 star 



4. Statistical properties of tlie sampie and 
phiysicai interpretation 

4. 1. Effective temperature scale 

To derive a reliable calibration for the effective temperature as a 
function of the spectral subclass for B dwarfs, we need enough 
stars for which we have accurate information for both parame- 
ters. Unfortunately, this prerequisite is not fulfilled for our sam- 
ple, as can be seen from the very large scatter around the ex- 
isting spectral-type-Teff-calibrations, presented in Fig.|9l The 
scatter is due to uncertain effective temperatures for some tar- 
gets and wrong spectral type designation for others. Indeed, for 
several stars, we encountered problems in the analysis, primar- 
ily because of the absence of Si lines for many fast-rotating late 
B type stars. These stars are too cool to have Si III present in 
their line spectrum, and the Si II lines are completely smoothed, 
so that they become hardly detectable. Even when they are still 
visible, the two lines of the Si II 41 28-41 30 doublet are so heav- 
ily blended that they cannot be used any longer. In these cases, 
the He I lines are the only lines that can be used to estimate the 
effective temperature, which are, therefore, not reliable enough 
to use in the derivation of an accurate temperature calibra- 
tion. On the other hand, several stars for which we were able 
to derive accurate values for the effective temperature turned 
out to have completely wrong spectral types. We caution the 
use of spectral types as they are quite often derived from low- 
resolution and low-quality spectra. Flags and remarks on the fit 
quality of individual targets were added in Table[3j 
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Figure 6. Example of the line profile fits for a 'cool' giant which has, meanwhile, been observed by CoRoT (MigUo et al., 
in preparation): HD 181440 (B9 III). Despite the low quality of the data, and the fact that we have no more Sim left in this 
temperature region, we are still able to obtain a satisfying solution. 



4.2. The stars in the (Tes, log gc)-cliagram 

Fig. [To] shows the (Teff, log g£.)-diagram of the analyzed 
GAUDl stars, along with the position of instability strips. 
The la tter were computed with the evolutionary code 
CLES ('Scuflaire etall l2008h and the pulsation code MAD 
(iDupr e t et al. 2002 ), respectively. These strips were computed 
for a core-ov ershoot parameter o f 0.2 times the local pressure 
scale height ( Miglio et al. 2007allbl) . It was indeed found from 
seismic modeling of B type sta rs that core o vershooting occurs, 
on average with such a value (lAertsll2008h . We thus compare 
the position of the GAUDI stars with the most recent models 
tuned by asteroseismology. 

For 40% of the GAUDI stars (26 out of 66), we have in- 
formation on the luminosity class. All but one (i.e. 9 out of 10) 
with known luminosity class II or III are indeed situated beyond 
the Terminal Age Main Sequence (TAMS hereafter), so we 



confirm them to be giants. The only exception is HD 49643, a 
star of spectral type B8 Illn, i.e. with spectroscopic hnes which 
may originate from nebulosity. The star is a fast rotator with a 
V sin / of almost 300 km s ' . The absence of Si lines due to the 
high rotation implies less reliable estimates for T^ff and \ogg 
than for the bulk of the sample stars. Thus, its derived gravity 
of 3.70 dex might be compatible with its previous classification 
as a giant. 

Three of the 14 stars with luminosity class IV or V turn 
out to be giants, i.e. they have - within the error bars - a sur- 
face gravity \oggc below 3.5. It concerns HD 43461 (B6V), 
HD 50251 (B8V), and HD 182198 (B9V), of which HD 43461 
is a fast rotator and has, therefore, less reliable parameters 
due to the absence of Si lines. Regarding the other two stars, 
there is a certain possibility that they are fast rotators as well, 
observed almost pole-on. Due to centrifugal effects, the stel- 
lar surface would become distorted, and effective tempera- 
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Figure?. Example of the line profile fits for the only hot supergiant in our sample: HD 52382 (Bllb). While the peak of 
Ha is reasonably well reproduced, the blue-ward absorption trough cannot be fitted at the available grid combinations of wind 
strength parameter and wind velocity field exponent. Note that no interpolation is performed for the se quantities, and that the 
wind is considered as unclumped, which might explain the obvious mismatch ('e.g.. lPuls et al. I l2006t in particular their Fig. 7). 
The Balmer lines show bumps in the wings. It is not clear whether this is real and due to the strong wind, or if this is a spectral 
artefact. It also arises in the only cool supergiant in our sample. The interpolated Si abundance (log n(Si)/n(H) - -4.67) is higher 
than the closest grid abundance (log n(Si)/n(H) - -4.79, plotted here), which explains the observed discrepancies in the Si line 
cores. In dotted lines, we show the line profiles as they would appear with the interpolated parameters. 
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Figure 8. Example of the line profile fits for an intermediate hot dwarf: HD 44700 (B3 V). The strong lines on the blue side of 
Si III 4716 and Sim 4819 are easily mistaken for the Si lines themselves. This is the reason why we always indicate the exact 
position of the transition during the full preparation process. The misfit of He I 4010 and 4140 is due to incomplete broadening 
functions of these lines. These lines were not used during the fitting procedure, but only as a double check afterwards (see 
Appendix lB . 11 1. We observe a similar behavior in several other stars. All other lines, even the weakest, fit very well. The shape 
of He I 447 1 is even perfect. 
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Figure 9. Effective temperature as a function of spectral type. 
We compare the position of our sample stars to existing tem- 
perature scales: for G alactic dwarfs, for spectral types up to B3 
( Trundle et al|2007 , solid blac k line), for dwarfs in the entire B 

, dashe d black line), for Galactic 
2007bi dotted 



type domain (Crowther 1998 
B supergiants CLefever et al 



Galactic and SMC B supergiants jMarkova & Pulsll2008 l solid 
grey line). Due to wind blanketing effects, the calibrations for 
the B supergiants lie at lower temperatures than the calibrations 
for B dwarfs. The supergiants (I.e. 1) are indicated by diamonds, 
the (sub)giants (I.e. 11 and 111) by asterisks and the (sub)dwarfs 
(I.e. IV and V) by squares. Filled circles indicate the stars for 
which we do not have any a priori information on luminosity 
class. It can be seen that spectral types taken from the literature 
can be quite inaccurate (see text). 



ture (due to gravity darkening) and surface gravity increase 
towards the pole. If observed pole-on, both quantities might 
be underestimated with respect to their average values. On the 
other hand, also the main sequence gets extended compared 
to standard models whe n rotation is taken into account (e.g., 
Maeder & Mevnetl2000l) . Simil ar results to ours were obtained 
by Hempel & Holwegen ( 2003 ). who found 6 out of their 27 
sample stars with spectral type from B6V to B9V to have a 
log § below 3.5. 

Depending on the input physics, stellar evolution theory 
shows that roughly 10 - 20% of the B8-B9 star s are in the gi- 
ant phase, but only 2 to 5% of B0-B3 stars (e.g., |Prialnikll2000t 
Scuflaire et al.ll2008h . >From the 40 stars without luminosity 
class treated by us, 1 1 turn out to be beyond the main sequence. 
From the 14 dwarfs, 3 turn out to be giants as well, and from the 
10 giants, one turned out to be a dwarf instead. Thus our sam- 
ple contains 23 giants, 19 of which have spectral type B8-B9. 
As our sample contains in total 51 stars of spectral type B8-B9, 
this means that the percentage of late type giants in our sample 
(37%) is about twice the expected one. Only five stars in our 
sample have spectral types earlier than B3, one of which is a 
giant, which is too low number statistics to verify the expected 
percentage. We found 20 of the 40 stars without luminosity 
class to have vsin/ above 90 kms and for 8 of them even 
V sin / > 150 km s ' . This fast rotation may again explain why 
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Figure 10. (reff, \og gc) diagram of the analyzed GAUDl B 
stars, with the corresponding 2cr-error bars. Symbols are the 
same as in Fig.|9] The dotted line represents the ZAMS, and five 
initial ZAMS masses - in Mq - have been indicated. The most 
recent theoretical instabiUty domains for the /3 Cephei (thick 
solid line) and the SPB stars (dashed lines) for main-sequence 
models, with a core overshooting value of 0.2 times t he lo- 
cal pressure scale height, are shown ( Miglio et alj|2007allbl) . to- 
gether with the instability domains for post-TAMS models with 
{ = 1 (grey dotted) a nd { = 2 (black dotted) g-modes computed 
bv lSaio et al. I (l2006h . The TAMS is shown as the low-gravity 
edges of the SPB and /3 Cephei instability domains around 
log gc = 3.5. 



we observe twice the number of expected giants when com- 
paring with standard models. Rotation was also found to be a 
necessary ingredient in evolutionary models in order to explain 
the number of observed giants in t he well-studied open clusters 
h and ;t'Per jVrancken et al ]|2000[ in particular their Fig. 3). 



4.3. The macroturbulent and rotational velocities 

When comparing the rotational broadening with macroturbu- 
lence, we observe a systematic trend of increasing Vmacro with 
increasing v sin / for those stars for which some pulsational be- 
havior may be expected, i.e., for the stars where Vmacio cannot 
be neglected (filled circ les in Fig.fTTTi. This is consistent with 
Markova & Pulsl (120081 and references therein), who found 
that, in almost all cases, the size of the macroturbulent velocity 
was similar to the size of the rotational velocity. On the other 
hand, they also found vsin/ and Vmacro to decrease towards 
later subtypes, being about a factor of two lower at B9 than at 
BO. 5, and that, in none of their sample stars, rotation alone was 
able to reproduce the observed line profiles. We cannot confirm 
either of both statements, as we have found, on the contrary, 
many cases with zero macroturbulence, and many late B type 
stars with very high macroturbulent velocities. The trend con- 
necting V sin / and Vmacio suggests that both broadening mech- 
anisms are difficult to disentangle. All Vmacro' -values seem to 
be at least one third of vsin ; , as can be derived from the posi- 
tion of the filled circles above the dashed line in Fig.[TT] From 
a certain projected rotational velocity on (i.e., around vsin; 
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Figure 11. Comparison between the macroturbulence and 
the projected rotational velocity v sin / . The solid line is the 
1-1 relation, while the dashed line represents the relation 
Vmacro/vsin; = 1/3. Those stars which we consider as non- 
pulsating (i.e., which have a negligible Vmacro) are indicated 
by open circles. The filled circles represent suspected pulsators 
(i.e., they have a significant Vmacro)- 



= 120 kms dotted vertical line), it is impossible to separate 
both effects, as the pulsational behavior (represented by v,„acio) 
completely disappears in the projected rotational broadening. 
Following the earlier argumentation that stars which have a 
non-negligible macroturbulent velocity may have the largest 
pulsational amplitudes, we investigated the behavior of Vmacio 
in the (Tes, logg)-plane (see Fig.fTSIi. On a global scale, the 
Vmacro-values show a random spread. The position of some in- 
dividual cases is intriguing, e.g., the supergiant HD 52382 (log 
Teff ~ 4.36, \ogg ^2.71). The star's position is compatible with 
the occurrence of gravity modes excited by the opacity mech- 
anism in evolved stars and may bel ong to a recently fo und 
class of pulsati ng B-type supergia nts ( Lefever et al.l2007bl) . As 



pointed out by lAerts et al.l (l2009h . these gravity-mode oscilla- 
tions may result in pulsational line broadening, hence the high 
Vmacro-value wc find. HD 48807 (B7 lab) and HD 48434 (BO 111) 
also show macroturbulence and may belong to the same class. 
Some of the dwarfs in the sample show macroturbulence, e.g., 
HD 48215 (B5 V), HD 177880 (B5 V) and HD 50251 (B8V) 
might be g-mode pulsators in the SPB instability domain. Many 
stars with \ogg below 3.5 show a reasonable to large vmacro, 
which may imply them to be pulsators. It is noteworthy that all 
instability computations for mid to late B stars so far have been 
computed for non-rotating models and were stopped artificially 
at the TAMS from the argument that no star is expected to be 
found in the Hertzsprung gap. B8 to B9 stars, however, cross 
this gap at a much lower pace than stars with spectral type ear- 
lier than B7 and our observational spectroscopic results point 
out that it would be worthwhile to extend the instability strip 
computations past the TAMS. 

Finally, the magnitude of the microturbulence seems to be 
randomly distributed in the (Teff, log g)-diagram, and in no way 
related to the effective temperature nor the surface gravity (see 
Fig.inil. 
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Figure 12. Same as Fig.[TOl but representing the values of 
i^macio for the GAUDl B star sample. The size of the symbols is 
proportional to the Vmacro-value. 
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Figure 13. Same as Fig. [121 but the size of the symbols is now 
proportional to the microturbulent velocity ^ of the GAUDI B 
star sample. 



5. Future use of our results 

The GAUDl database was set up during the preparation of the 
CoRoT space mission. We applied our developed automated 
tool for the determination of the fundamental parameters to nu- 
merous B stars in this database. A recent independent study 
used the same database to determine the abundances of the 
89 B6-B9 stars in GAUDl, using a different approach based 
on LTE atmosphere models without wind and LTE occupation 
numbers. We refer to Niemczura et al. ( 20091) for the underly- 
ing assumptions, but mention here that a mild depletion of both 
Fe and Si, with average values of the iron and silicon abun- 
dances of 7.13+0.29 dex and 7.22+0.31 dex, respectively, was 
found for their sample. This result was obtained from spectrum 
synthesis in the LTE approximation, by fixing the effective tem- 
peratures from photometric calibrations, the microturbulence at 
2kms"' and ignoring any macroturbulence in the line profile 
fits. The mild Si depletion is in agreement with our results for 
the Si abundance of the late B stars treated by us. 
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Meanwhile the first data of CoRoT were released to the 
public. Two of our sample stars were already observed by 
the satelhte. It concerns HD 181440 (B9 111) and HD 182198 
(B9 V). The spectra of these stars were available in the GAUDI 
database and their parameters are listed in Table[3] Their seis- 
mic analysis is under way (Miglio et al., in preparation). 

We also point out that one of the primary targets of CoRoT 
is HD 180642 (B1.5 II-III), which is a mildly N-enriche d, ap- 
parent l y slowly rotating 13 C ephei star, discovered by lAerts 
(l2000l) . iMorel & AertsI (l2007h determined its fundamental pa- 
rameters from high-resolution spectroscopy and derived T^s = 
24500+ 1000 K and logg = 3.45 + 0.15,^= 12 + 3kms-i.By 
means of AnalyseBstar, we found the model with Teff = 24000 
+ 1 000 K, logg = 3.40 + 0.10 and ^ = 13 ± 2 kms'^o be the 
best fitting one, so our results agree very well with theirs. We 
thus classify this star as a pulsator very close to the TAMS. It 
will be very interesting to know if this spectroscopic result will 
be confirmed by the seismic modeling of oscillation frequen- 
cies derived from the CoRoT space photometry, because this 
would imply one of the few /3 Cephei stars very close to the end 
of the core-hydrogen burning stage. 

At the time of launch, the CoRoT space mission was nom- 
inally approved to operate through 2009, but, given its excel- 
lent performance, the operations will probably be prolonged. 
It is therefore to be expected that several more B-type stars in 
GAUDI will be observed by the mission in the near future. 
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Appendix A: Formal convergence tests for 
synthetic FASTWIND spectra - some details 

Before applying AnalyseBstar to real, observed spectra, we 
first tested whether the method is able to recover the input 
parameters of synthetic spectra. For this purpose, we have cre- 
ated several synthetic datasets in various regions of parameter 
space. We will not dwell on discussing all of them here, but we 
have chosen three specific examples, each representative for a 
different type of star: 

- Dataset A: a BO. 5 I star with a rather dense stellar wind and a 
strongly enhanced helium abundance, 

- Dataset B: a B3 III star with a weak stellar wind and a 
depleted Si abundance 

- Dataset C: a B8 V with a very thin stellar wind. 

Each dataset has been convolved with both a rotational and 
a macroturbulent broadening profile. The projected rotational 
velocity adopted is 50 kms"' for each dataset, characteristic 
for a 'typical' slow rotator. In the case of low and interme- 
diate resolution spectra, there is also the possibility to carry 
out an additional Gaussian, instrumental convolution. This is, 
however, not necessary for spectra with such high resolution 
as FEROS and ELODIE. Artificial (normally distributed) noise 
was added to mimic a real spectrum with a mean local SNR of 
150, which is more or less the minimal local SNR obtained for 
the GAUDI sample. 

Finally, we have gone through the full process of the prepa- 
ration of the spectra as if it were real data, i.e., defining the EW 
of all available lines, measuring the SNR to account for the 
errors on the EW and fixing the projected rotational velocity. 
Since we are dealing with synthetic data, obviously no normal- 
ization was required. 

Table lA. Il lists the input parameters for the three synthetic 
datasets as well as the output parameters, obtained from the ap- 
plication of AnalyseBstar. For each dataset, the second column 
gives the derived 'interpolated' values, while the third column 
lists the parameters of the closest (best fitting) grid model. In 
the ideal case, both should be exactly the same as the input pa- 
rameters. This depends, however, on how well the equivalent 
widths were measured and minor deviations occur as expected. 
Besides these three cases, in which we started from a synthetic 
model which is one of the grid points, we have additionally cre- 
ated three synthetic datasets for which the parameters lie in be- 
tween different grid points. In analogy to datasets A to C, their 
profiles have been convolved with both a rotational and macro- 
turbulent broadening profile and artificial noise was added. The 
input parameters and the results of the analysis have been added 
to Table lA. ri as datasets D, E and F. 

A third set of models (G, H, I) consists of models computed for 
parameters not included in the grid, but within the grid limits, 
and this time without noise, with the aim to test the predictive 
power of our method, independent of the source of errors in- 
troduced by the noise. The rotational velocity applied for these 
test models is 30 km s and the macroturbulent velocity is 15 
kms"^'. We find an overall very good agreement of the closest 
grid model and the interpolated values with the input parame- 
ters. 



Although the models given in Table lA. ll were selected as 
the best fitting models, the program additionally came up with 
some other models, which agree with the input model within 
the errors, introduced by the artificial noise and the errors in 
the determination of the equivalent widths (datasets A to F) 
and/or the various interpolations (datasets D to I). For instance, 
for dataset I, AnalyseBstar came up with a second possibility 
with slightly different parameters, even closer to the input data: 
log^ = 1.7, ^ = 6.5 kms-', v,nacro = 14.6 kms-', N(He)/N(H) 
= 0.12, but with log n(Si)/n(H) = -4.49 (the other parameters 
are the same as the ones for the model in the table). As the Si 
abundance is overestimated compared to the input value, the 
Si lines are overall a bit too strong and the final fit is slightly 
worse than the one for the model in the table. 

In all nine cases, the input parameters are well recovered. 
Typically 10 to 150 models are selected and considered during 
the analysis cycle. The number depends not only on the choice 
of the initial parameters, but also on the accuracy of the mea- 
sured quantities and the applied procedure (hence, the temper- 
ature domain, see below). We considered numerous other test 
cases, which are not included in this text, but for which the in- 
put parameters were equally well recovered. Slight deviations 
from the input parameters are as expected. Also in the cool tem- 
perature domain (datasets C, F and I), where our alternative 
optimization 'method 2' is required (see Section lB3] l. we are 
still able to deduce reliable parameter estimates, at least for the 
synthetic models. The results for the Si abundance in datasets F 
clearly show the problem for estimating the abundances when 
only one ion is present: the estimation of the Si abundance was, 
for this particular test case, not equal to the input value. When 
we have only Si II lines, the results cannot be as secure as when 
we have two ionisation stages of Si. This problem has been ac- 
counted for by adopting somewhat larger errors on the derived 
parameters, for objects with Teff < 15 000 K (see Table[3]i. 
When the observed object is close to one of the borders of the 
grid, the method takes 'the border' as the closest match. This 
will, consequently, lead to larger errors in some of the other 
parameters. 

The derived v sin / -value, which is difficult to disentangle 
from the macroturbulent velocity, agrees very well with the in- 
serted values, irrespective of v,naci-o- This shows that the Fourier 
Transform method indeed allows to separate both effects. All 
together, this gives us confidence that our method is working 
reliably and that the procedure will also be able to recover the 
true physical parameters from real spectra. 

The macroturbulent velocity is the only fit parameter for 
which significant deviations arise, if Vmacro is low (< 30 
kms '). Larger macroturbulences are well recovered, how- 
ever Especially for datasets C, E and I, where Vmacro is only 
10-15 kms the deviation is large (of the order of 15 to 
30 kms '). To understand this discrepancy, we first had a look 
at the line profiles, which show that the fit is almost perfect for 
this Vjnacro (seeFig lATTT i. After verifying the fit quality, we also 
verified that the discrepancies do not arise from our applied 
procedure, by performing several tests on simulated spectra, in 
which we left Vmacro as the only free parameter. In all cases, the 
macroturbulent velocity was recovered very well, with devia- 
tions within 5 km s ' . The (sometimes significant) deviations in 
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^macro Can be understood as a compensation for differences be- 
tween the interpolated values of the 'observ ed' line profiles an d 
the closest grid model. Indeed, as shown bv lAerts et al. I (l2009h . 
small deviations in the wings of the Si profiles can result in 
wrong estimates of Vmacro and v sin / , because these two veloc- 
ity fields are hard to disentangle whenever one of them is small. 
Thus, the values derived for Vmaao should be treated with cau- 
tion. This does not affect the derivation of the other parameters, 
however, since a convolution with the Vmacio-profile preserves 
the EW of the lines. It is the EW, and not the line profile shape, 
which is used for the derivation of Teff, ^ and log n(Si)/n(H). 
The other parameters (gravity, mass-loss rate, and wind veloc- 
ity law) remain unaffected because they are deduced from the 
Balmer lines, which are not very sensitive to the Vmacro values 
due to the dominance of the Stark broadening. Thus we find 
a very good agreement between their input and output values 
(TablelAJ}. 



Appendix B: Methodology of AnalyseBstar 

In what follows, we give a detailed description of AnalyseBstar 
In Fig. IB. II a flowchart of the full program is drawn schemat- 
ically. The reader is advised to follow this context diagram 
throughout the further description of AnalyseBstar. 



B. 1. Preparation of the input 

No automatic method can fully replace the 'by-eye' procedure, 
and this is also true for AnalyseBstar. A few steps still require 
human intervention. This especially concerns the preparation 
of the spectra and the input for the main program, as will be 
shown in the following. 



Observation of peculiarities. The first thing to be done when 
starting the analysis is to inspect the spectrum regarding any 
'abnormalities', such as line behavior resulting from nebulae, 
disks, binarity or other peculiarities. Special features can and 
should be detected through eye inspection. 

Merging and normalization of the spectra. When dealing 
with echelle spectra, like the ELODIE and FEROS spectra in 
the GAUDI database, we first need to merge the spectral orders 
before the spectrum can be normalized. This is a very delicate 
and non-trivial process. The edges of the orders are often noisy, 
due to the decrease in sensitivity of the CCD near its edges, and 
are therefore removed. Spectral lines falling inconveniently on 
the edge of one spectral order or in the overlap region of two 
consecutive orders, thus resulting in a cut-off within the spec- 
tral line, are not reliable to derive accurate stellar information. 
They can, in the best case, only be used for a consistency check 
afterwards. After the spectra have been merged, they are nor- 
malized through interactive interpolation of continuum regions. 
This is done over a large wavelength coverage, where sufficient 
continuum regions are present, in order to avoid line corrup- 
tion which will propagate in the derivation of the fundamen- 
tal parameters. The normalization especially proves to be very 



important for the determination of the gravity, which depends 
sensitively on the Balmer line wings. 

Radial velocity correction. The spectra need to be corrected 
for radial velocity shifts. Shifting the spectrum to rest wave- 
length automatically has not been included in the current ver- 
sion yet. 

Line selection. As a standard procedure, we compute a num- 
ber of line profiles in the optical from the underlying model 
atmosphere. These standard line profiles can be used directly. 
Computing additional line profiles is, however, straightfor- 
ward. The selection of lines to be used during the spectral anal- 
ysis in AnalyseBstar is gathered from a file which contains a 
default (but easily editable) list of currently used optical lines 
(see below). We distinguish between optical lines that will be 
used for the prediction of the physical parameters during the 
spectral line fitting procedure and those that will not be used 
explicitly for reasons of uncertainties, either in the theoreti- 
cal predictions or in the observed spectrum (due to problems 
in merging or normalization). The goodness-of-fit of the lat- 
ter lines is only checked to investigate systematic differences 
between theory and observations, which can be useful input 
for considerations regarding further improvements of atomic 
data and/or atmospheric models. They are He, He 1 4010, 4120, 
4140 and Si IV 4950. The lines that will be used in the fit- 
ting procedure for their predictive power are the following: the 
Balmer lines Ha, H/J, Hy and H<5, the neutral and singly ionized 
He lines He I 4026, 4387, 4471, 4713, 4922, 6678, He II 4541 
and He II 4686, and the Si lines in their three different ioniza- 
tion stages Sill 4128-4130, Sill 5041-5056, Sim 4552-4567- 
4574, 4716, 4813-4819-4829, 5739, Si IV 41 16 and Si IV 4212. 
He II 4200 was excluded as it can only safely be used for O- 
type stars and not for B-type stars. At least for dwarfs, it only 
appears at temperatures above about 25 000 K, where it is still 
overruled by N II, while at hotter temperatures. He II 4200 is 
for 50% blended by Nlll^ 

Signal-to-noise measurement. To estimate the error on the 
equivalent width as accurate as possible, we need to determine 
the signal-to-noise ratio, SNR. Since the SNR can vary a lot 
with wavelength, we chose to calculate the local SNR for each 
line separately, by using a continuum region close to the line. 
This is done manually through the indication of the left and 
right edges of the continuum interval. The SNR is then given 
by the mean flux divided by the standard deviation a-{Fc) of 
the flux in this interval, i.e., 

SNR = Tjcr(F,). 



Observed equivalent widths and their errors (cf. Fig.| 
The equivalent widths of most He and Si lines can be measured 
by a Gaussian non-linear least squares fit to the ob served line 
profiles, using the Leve nberg-Marquardt algorithm ( Levenber3 
1944l:lMarquardtll963h . We thoroughly tested whether this was 



See |http://www.lsw.uni-heidelberg.de/cgi-biii/websynspec.cgi| 
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Table A.l. Input parameters for the synthetic models (IN) are compared to the actual output parameters (OUT) obtained through 
the application of AnalyseBstar to these synthetic data: the effective temperature, T^f^, the surface gravity, log g, the wind strength 
parameters, log Q (see below), the wind velocity exponent, fi, the He abundance, n(He)/n(H), the Si abundance, log n(Si)/n(H), 
the microturbulent velocity, ^, the projected rotational velocity, vsin/ , the macroturbulent velocity, Vmacro, the time elapsed in 
seconds during the run of AnalyseBstar, and the number of different models that have been checked throughout the procedure. 
The wind strength parameter log Q is given by log M/(Vco R*)'^, where M is the mass-loss rate, Voo the terminal wind velocity, 
and Rt the stellar radius. In our grid, we allow for 7 different values, represented by a character: log Q = -14.30 (O), -14.00 (a), 
-13.80 (A), -13.60 (b), -13.40 (B), -13.15 (C), -12.70 (D). 

Datasets A, B, and C are models which lie exactly on a grid point, whereas dataset D, E, and F are models which lie in between 
the gridpoints. Artificial noise was added to the synthetic line profiles in these two datasets, to simulate a real spectrum. The 
models for datasets G, H, and I also lie in between the gridpoints, but no artificial noise was added and a lower vsin i (30 km s ' 
instead of 50 km s"') was used in order to test the predictive capabilities of the method. The 'interpolated' values derived for 
the parameters, as well as the best fitting (i.e. closest) grid model are given. The parameters in italics are those for which no real 
interpolation was made, but for which the most representative grid value is chosen. The longer computation times for the cooler 
models are due to the different method applied in this range. 
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2344 






checked models 


10 






21 






128 








Dataset D 


Dataset E 


Dataset F 


Teff (K) 


21,120 


20,985 


21,000 


15,100 


15,030 


15,000 


11,880 


11,500 


11,500 


log,?(cgs) 


3.98 


4.0 


4.0 


1.83 


1.80 


1.8 


2.43 


2.3 


2.3 


log Q (char) 


b 


b 


b 


A 


A 


A 


a 





O 


P 


1.42 


2.0 


1.2 


2.8 


3.0 


3.0 


1.02 


0.9 


0.9 


n(He)/n(H) 


0.14 


0.11 


0.10 


0.08 


0.09 


0.10 


0.18 


0.20 


0.20 


log n(Si)/n(H) 


-4.85 


-4.80 


-4.79 


-4.23 


-4.16 


-4.19 


-4.49 


-4.19 


-4.19 


f (kms-') 


12 


15.7 


15 


11 


10.6 


10 


9 


6.3 


6 


vsin/ (kms"') 


50 


52 ±6 


52 ± 6 


50 


50 + 7 


50 + 7 


50 


54 ± 1 


54 ± 1 


Vraacro (km S"') 


40 


33 ±5 


33 ±5 


10 


39 ±7 


39 ±7 


70 


67 ±7 


67 ±7 


time (s) 


174 






392 






1254 






checked models 


10 






26 






95 








Dataset G 


Dataset H 


Dataset I 


Teff (K) 


24,800 


24,650 


25,000 


18,600 


18,392 


18,500 


10,910 


11,000 


11,000 


logg(cgs) 


3.33 


3.3 


3.3 


2.78 


2.8 


2.8 


1.72 


1.8 


1.8 


log Q (char) 


B 


B 


B 


O 


A 


A 


C 


C 


C 


P 


1.1 


1.2 


1.2 


1.0 


0.9 


0.9 


2.10 


2.0 


2.0 


n(He)/n(H) 


0.11 


0.12 


0.10 


0.14 


0.13 


0.15 


0.11 


0.09 


0.10 


log n(Si)/n(H) 


-4.52 


-4.48 


-4.49 


-4.25 


-4.22 


-4.19 


-4.67 


-4.79 


-4.79 


f (kms"') 


13 


13.5 


12 


4 


3.3 


3 


7 


8.6 


10 


vsin; (kms~') 


30 


30 ± 1 


30 ± 1 


30 


30 ± 1 


30 ± 1 


30 


30 ± 1 


30 ± 1 


Vraacro (km S"') 


15 


12 ±4 


12 ±4 


15 


20 ±7 


20 ±7 


15 








time (s) 


2119 






2284 






356 






checked models 


66 






67 






15 
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HEI4- 1 20 




Figure A.l. Example of the fit quality for dataset E. The shght discrepancy in the cores of the Si lines arises from the difference 
in Si abundance between the interpolated value (log n(Si)/n(H) = -4.16) and the closest grid value (log n(Si)/n(H) = -4.19) 



a justifiable approach by comparing, for multiple line profiles 
of stars with completely different properties, the equivalent 
widths computed from a Gaussian fit to the observed line pro- 
files with those computed from a proper integration of the line 
pixel values. This comparison resulted in a minimal difference 
between both approaches. Some lines obviously cannot be fit- 
ted by a Gaussian, e.g., the He fines with a strong forbidden 
component in their blue wing or Stark-broadened Unes such 
as those from hydrogen. The EW of such lines are determined 
by proper integration. In all other cases, where the line profile 
is well-represented by a Gaussian fit, the the observed EW is 
given by the integral of the Gaussian profile. 
The main source of uncertainty in the EW determination is in- 
troduced by the noise and its influence on the exact position of 
the continuum within the noise. To account for this we multi- 
phed the continuum by a factor (1 ± 1/SNR), performed a new 
Gaussian fit through both shifted profiles and measured the dif- 
ference in EW. Note that, by considering a constant factor over 
the full line profiles (using the signal at continuum level), we 
implicitly give an equal weight to each wavelength point, even 



though, in reality, 1/SNR may be slightly larger in the cores, 
since it is proportional to VI IS , with S the obtained signal at 
these wavelengths. When the user deems it necessary, (s)he can 
apply a correction for the continuum level by a local rerectifi- 
cation, e.g., in the case that the normalization performed over 
a large wavelength range is locally a bit offset. If the offset 
is different on the blue and red side of the line, then the fac- 
tor should be wavelength dependent. Therefore, we allow for 
a linear rerectification by a factor aA + b. Another source of 
uncertainty in the EW is introduced by using a Gaussian fit to 
estimate the EW. Each of the above factors have been included 
in the total error budget of the EW. 



Determination of vsin/ (see Fig. lB.3b . The determination 
of the projected rotational velocity is the last step in the 
preparation of the input for AnalyseBstar. We use t he sem i- 
automatic tool developed by ISimon-Diaz & Herrerol (l2007b to 
derive the projected rotational velocities (see Fig. lB.3l l. We re- 
fer to their paper for a thorough discussion. We use the fol- 
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Section B.8 




SNR, EW ± AEW, vsini 
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STARTING VALUES 
Section B.2 




Si EW line ratios: 
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Section B.3 (1) 



EW of Si and He: 
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Section B.3 (2) 
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Section B.3 (1) 



Fit of Si line wings: vmacro 

Section B.4 
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Fit of Balmer lines: logg 

Section B.5 



Fit of Ha profile: log Q, B 

Section B.6 
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Figure B.l. Context diagram of AnalyseBstar. For a detailed explanation, see text. 



lowing, least blended, metallic lines for the determination of 
vsin/ : Sill 5041, Si III 4567, 4574, 4813, 4819, 4829, 5739, 
Cii 4267, 6578, 6582, 5133, 5145, 5151 and On 4452. We 
carefully checked visually that blended lines were not taken 
into account. The projected rotational velocity v sin ; of the star 
is calculated as the mean of the values derived from each indi- 
vidual line, and its uncertainty as the standard deviation. 

B.2. Fit parameters and their starting values 

Once all preparation is finished, the automatic procedure to 
obtain accurate values for all physical fit parameters can be 
started. This procedure follows an iterative scheme as outlined 
in Fig. IB. II In each iteration, the fundamental parameters are 
improved in the following order. First of all, the effects of the 
effective temperature, the Si abundance and the microturbu- 
lence on the Si lines are separated. Then, the He abundance 
is fixed using all available He lines. Next, the macroturbulent 
velocity is determined from well-chosen, user supplied Si lines. 



Then, the surface gravity is determined from the wings of Hy, 
W5 and, in the case of weak winds, also H/?, after which the 
wind strength log Q (Q = A^f/(vco ^»)'^, with M, the mass loss, 
Voo, the terminal wind velocity, and the stellar radius) and 
the wind velocity exponent yS are determined in parallel using 
Hq'. Obviously, to start the iterative procedure, we need an ini- 
tial guess for each of these parameters. This initial value can 
either be user supplied or standard. If no value is set by the 
user, then the following initial values are considered. The ini- 
tial effective temperature will be determined from the spectral 
type of the star Immediately after the rotational and macro- 
turbulent velocities have been determined, a start value for the 
surface gravity is derived by running the subprocedure for the 
determination of the gravity (see Section lB3T l. For the He and 
the Si abundance, we have taken the lowest value as start value. 
We initialize the wind parameters at the lowest values, which 
means negligible wind, at logQ = -14.30 and fi = 0.9. These 
are reasonable assumptions for dwarfs, which comprised the 
largest part of our sample. These starting values can easily be 
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3 f^T-VVfi »-5~s 



402^ 4026 4028 4030 



Figure B.2. Illustration of the EW determination for He 1 4026. 
After the manual identification of the wavelength interval of the 
spectral line (thick black profile), a Gaussian fit to the observed 
line profile is made to determine the EW of the line (thick grey 
profile). Also indicated are: the center of the line (vertical line) 
and the 'shift' in flux, upwards and downwards (dashed line 
profiles), used to account for the noise level in the determi- 
nation of the EW. Also to these line profiles a Gaussian fit is 
made, but these are omitted here for the sake of clarity. 




0,0 



0.10 1.00 
sigma(cycles/A) 



Figure B.3. Illustration of the determination of v sin / from the 
first minimum in the Fourier transform of the selected li ne pro- 
file, as implemented and described by iSimon-Diaz & Herrerol 
(l2007h . The user can make several attempts and finally decide 
which V sin / gives the best match. The difference in the slope 
of the first decay gives an indication of the macroturbulence. 
The Fourier transform of the observed line profile is indicated 
in black. The user indicates the first minimum in the black pro- 
file, which gives the projected rotational velocity. The Fourier 
transform of the rotational profile at this vsini is indicated in 
grey. 



adapted when samples of stars with different stellar proper- 
ties are aimed at. The macroturbulent velocity is initially zero, 
whereas for the microturbulent velocity we took a medium 
value of all considered possibilities, i.e., 10 kms"'. The con- 
vergence speed depends on how far away the initial values are 
from the final solution. 



1 .0 
0.5 
0.0 
-0.5 
-1.0 



1 .5 
1 .0 
0.5 
0.0 
-0.5 
-1.0 



• f - 



Tetf (K) 



• I 



Teff (K) 

Figure B.4. Synthetic simulation of the effect of the micro- 
turbulence (top) and the Si abundance (bottom) on the log- 
arithmic EW ratios of Si lines of different ionization stages. 
This example gives the (synthetic) EW ratio of Si III 5739 to 
Sill 4130 against a selected range of effective temperatures, 
for a fixed surface gravity (log g = 2.5) and wind parameters {M 
= 0.144.10"^Mo/yr, = 500 kms"', /? = 0.9). The ratios are 
evaluated for the different possibilities of the microturbulence, 
for a Si abundance fixed at -4.79 (upper panel; ^=3,6, 10, 12, 
15 and 20 km s"' ; triangle up, square, diamond, circle, asterisk, 
triangle down, respectively) and for the different possibilities 
of the Si abundances, for a microturbulence fixed at 10 kms"' 
(log n(Si)/n(H) = -4.79, -4.49 and -4.19: circle, triangle up and 
square, respectively). The horizontal lines show the observed 
EW ratio and its corresponding uncertainty region. Acceptable 
EW ratios, which fall within these boundaries, are indicated in 
grey. 



B.3. Determination of the effective temperature, 
microturbulence and abundances 

The Si lines serve multiple purposes. By using a well-defined 
scheme (somew hat similar to the conventional 'curve of 
growt h' method ( Gray 19761) and described, e.g., by Urbanejal 
20041) . we are able to separate the effects of the effective tem- 
perature, the Si abundance and the microturbulence on the line 
profiles, and derive an accurate value for them. Our method 
automatically determines from the observed spectrum which 
lines will be used for this purpose. The lines should be well 
visible (i.e. clearly distinguishable from the noise in the con- 
tinuum) and the relative error on their equivalent width is typ- 
ically below 10 to 15%, depending on the combination of the 
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temperature and the specific line considered. Whenever the rel- 
ative error on its equivalent width exceeds 75%, the spectral 
line will be excluded in the following procedures. Such large 
errors are a rare exception and only occur for the weakest lines, 
which almost disappear into the noise level. 
Basically, we can discern two different cases, namely when 
multiple and consecutive ionization stages of Si are available, 
or when there is only one ionization stage of Si. Each needs a 
different approach, as we will explain now. 

Method 1 : In case there are multiple ionization stages of 
Si available, we can, to first order, put aside the effect of the 
Si abundance by considering the EW ratios of two different 
Si ionization stages, as the Si abundance affects all Si lines in 
the same direction. For each effective temperature point in the 
grid, the ratios of the observed equivalent widths of Si IV to 
Si III and/or Si III to Si II are compared to those obtained from 
the model grid, given a certain log g and wind parameters (see 
Fig. lB.4b . For each combination of Si iv/lll and/or Si Ill/ll lines, 
this results in a range of 'acceptable' effective temperatures 
for which the observed EW ratio is reproduced, within the ob- 
served errors (indicated as grey symbols in Fig. IB. 4b . From the 
different line ratios, slightly different temperatures may arise. 
Under the assumption that log g and the wind parameters are 
perfectly known, the combination of all these possibilities con- 
stitutes the set of acceptable effective temperatures. 

Once we have a list of acceptable temperatures, we can de- 
rive for each temperature the best microturbulent velocity and 
the best Si abundance, by analyzing the variation in equivalent 
width of different lines as a function of these quantities. In what 
follows, we describe the full process, step by step; 

- Step 1: Deriving the abundances for each microturbulence, 
line by line 

In this first step, the observed EW of each Si line is com- 
pared to the set of theoretically predicted EWs correspond- 
ing to the available combinations of Si abundance and mi- 
croturbulence in the grid (see Fig. IB. 5b . For each microtur- 
bulent velocity (and for each Si line), we look for the range 
of Si abundances that reproduce the observed EW, within 
the observed errors. This is done through linear interpola- 
tion. Note that this gives us only an approximative estimate 
of the Si abundance. Indeed, as can be seen in Fig. IB. 51 the 
change in EW with increasing microturbulence is marginal 
for weak lines, but increases for stronger spectral lines, in- 
dicating that we may no longer be in the linear part of the 
curve of growth. 

- Step 2: Deriving the mean abundance for each microturbu- 
lence and the slope of the best fit 

In a second step we compare, for each microturbulence, the 
abundances (with the derived uncertainties) found in step 1 
to the observed EW of each line (see Fig. IB. 6b . Since a star 
can only have one Si abundance, we should find the same 
Si abundance from all the different Si lines, i.e. the slope of 
the best fit to all the lines should be zero. From the differ- 
ent abundances derived from each line, we can calculate the 
mean abundance for each microturbulence. 

- Step 3: Deriving the 'interpolated' microturbulent velocity 



In a third step, we investigate the change of the slope when 
varying the microturbulence (see upper panel Fig. IB. 7b . 
Through linear interpolation, we find the microturbulence 
for which the slope would be zero (i.e. for which the Si 
abundance derived from each line separately would be the 
same). This is the estimated 'interpolated' microturbulence. 

- Step 4: Deriving the 'interpolated' abundance 

Using now the relation between the microturbulence and the 
mean abundance derived in step 2, we interpolate to find 
the estimated 'interpolated' Si abundance at the estimated 
'interpolated' microturbulence found in step 3 (see lower 
panel Fig. IB. 7b . 

- Step 5: Deriving the 'interpolated' effective temperature 
Now that we obtained the estimated 'interpolated' values 
for the microturbulence and Si abundance, we can go back 
to Fig. IB. 41 and interpolate in two dimensions to derive a 
value for the effective temperature, which reproduces the 
observed EW. The mean of the effective temperatures from 
each line ratio will then be accepted as a value for the 'in- 
terpolated' effective temperature of the star. 

- Step 6: Determination of the 'closest' grid values 

For the next steps in the automatic procedure, we will need 
the closest grid values to these estimated 'interpolated' val- 
ues rather than the real values themselves. These closest val- 
ues will constitute a new entry in the list of possibilities. If 
the real value for the Si abundance or microturbulence falls 
exactly between two grid points, both are added to the list. 

We repeat these steps for each 'possible' grid temperature, ini- 
tially derived from Fig. lB.4l (grev symbols), and obtain in this 
way a set of new possible solutions that optimally reproduce the 
Si lines. Temporarily fixing the values for T^^, ^, log n(Si)/n(H) 
in the way described above leads to a new estimate for the He 
abundance, denoted as n(He)/n(H). As for the Si abundance, we 
only consider three different values. We apply a similar method 
as for the determination of the Si abundance, in the sense that 
we derive the best-suited abundance from each line separately, 
and take the mean value as the 'interpolated' value. Figs lB.81 
and IB. 91 illustrate this process. Note that we intrinsically as- 
sume that the microturbulence is the same throughout the at- 
mosphere, i.e. that there is no radial stratification. In this way, 
the microturbulence, necessary to account for the broadening in 
the He lines, will be assumed to be the same as the one derived 
from the Si lines. Note also that the uncertainty in the derived 
He abundance will be larger when the He abundance is lower 
than 0.10. Indeed, in this case we are forced to extrapolate to a 
region where it is unclear how the dependence of the equivalent 
width with abundance will change. The decrease towards abun- 
dances lower than solar may be steeper or follow a logarithmic 
trend, in which case the predicted values would be underesti- 
mated. In this case we can only state that the solar abundance 
is an upper limit of the true abundance. Values lower than the 
primordial He abundance of ~0. 10 can no longer be considered 
physical, except when diffusion effects start to play a role and 
He settles. In that case, the He abundance observed at the stellar 
surface could be lower. Such chemically peculiar stars are re- 
ferred to as He weak stars and are usually high-gravity objects. 
Helium settling was not taken into account in our models. 
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Figure B.5. Step 1 (for each Teff and for each line): For each of the five considered microturbulent velocities (5 dash-dotted lines, 
from low EW to high EW: 3, 6, 10, 12, 15 kms '), we derive the Si abundance that reproduces the observed EW (horizontal 
grey lines) of the spectral line through linear interpolation. The Si abundances and the corresponding upper and lower limit are 
represented by the crosses. The displayed example is a synthetic simulation. Only a selection of Unes is shown. Usually many 
more lines can be used. In this example, no Si IV was "observed", so only Si II and Si III were considered. 
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Figure B.7. Steps 3 & 4: the 'interpolated' microturbulent ve- 
locity and the accompanying 'interpolated' Si abundance (in- 
dicated in grey) are derived from the position where the slope 
of the best fit (step 2) is zero. The Si abundances, given in the 
lower panel, are the mean abundances derived in step 2. The 
displayed example is a synthetic simulation. 



Method 2: In case there is only one ionization stage 
of Si available (mostly for late B-type stars, for which we 
only have Si II, but also for few hotter objects around Teff = 
23 000 K ± 3000 K with Si III) only, we can no longer use 



the ionization balance to derive information about the effec- 
tive temperature, independent of the abundances of He and 
Si, and independent of the microturbulence. Indeed, we have 
only two known factors (in the case of late B-type stars: EW 
of He I and EW of Si ll) to derive four unknown parameters 
(T^eff, ^, log n(Si)/n(H), n(He)/n(H)). Fortunately, we can con- 
strain the effective temperature quite well. For instance, in the 
cold B-type domain. He I is very sensi tive to changes in ef - 
fective temperature (see, e.g.. Fig. 2 in Lefever et al. 2007ah . 
which means that we can use the joint predictive power of He I 
and Si II to derive a set of plausible values for Teff. This is done 
by comparing the theoretical EWs of all combinations of Si 
abundance, microturbulence and effective temperature (within 
a range ± 4 000 K around the currently investigated Tgif ) with 
the observed EWs, and selecting the 'best' combinati ons from 
an appropriate log-likelihood function (for details, see iLefevei 
i2007i) . This method does not only work for the cool domain, 
but also for those hotter objects where only Si III is present, in 
this case because the EW of Si III varies much faster than the 
EWofHel. 

For each of the derived plausible Tgff values, we still have 
to find the corresponding Si (from the Si lines) and He abun- 
dances (from the He lines), and the microturbulence (from 
the Si and/or He lines). As the He lines are liable to Stark 
broadening, we use only Si lines to derive the microturbulent 
broadening. This means that we essentially have to derive only 
three unknown parameters from two known data: two from the 
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Figure B.6. Step 2 (for each microturbulence): we plot the Si abundances and their uncertainties for all 6 lines for which we 
derived these values in step 1, as a function of the observed EW. The least squares fit to the abundances are shown in grey. 
The microturbulent velocity for which the slope of this fit is zero (i.e. equal abundance from each line) gives an estimate of the 
'interpolated' microturbulent velocity. The displayed example is a synthetic simulation. 
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Figure B.8. For each He line, the most suited He abundance is derived through interpolation. The theoretically predicted values 
are shown as filled circles, the interpolated values as crosses. The horizontal lines indicate the observed equivalent width and the 
observed errors. The displayed example is a synthetic simulation. 



Si lines (log n(Si)/n(H) and ^), and one more from the He lines similar), which still leaves us with one free parameter for which 
(n(He)/n(H), assuming that the microturbulence in Si and He is we will have to make an assumption in order to be able to fix 
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the other two. We have chosen to 'fix' the Si abundance, as it 
has no direct influence on the He hnes. However, we do con- 
sider each of the three possibilities for the Si abundance in the 
grid, i.e., we consecutively consider the three cases in which 
Si is depleted, solar and enhanced. This allows us to determine 
also the microturbulent velocity, which is consequently used to 
fix the He abundance from the He lines. 

To finally decide which of the various combinations 
matches the observations best, a twofold check is applied, (i) 
all profile-sets are inspected by eye (this, again, compromises a 
fully automatic procedure), since parameter combinations lead- 
ing to obvious mis-fits can be clearly excluded on this basis, 
(ii) Again, a log-likelihood procedure is used to find the best 
matching model among the list of final solutions, combining 
the likelihoods of the EWs of the Si /He lines and o f the profile 
shapes of all lines, respectively (see Lefeverll2007 ). 



Even though, with this procedure, we are able to derive rea- 
sonably good estimates for most of the stellar/wind parameters, 
we have no means to restrict the Si abundance in a reliable way 
as it is possible within 'method 1'. The outcome of maximizing 
the log-likelihood can give us only a general idea of the abun- 
dance, but still leaves us with a larger uncertainty, which has 
been accounted for in the adopted error budget (see Table|3]i. 




Figure B.9. The finally provided He abundance is calculated 
as the mean over all lines. The abundance derived from each 
line separately is represented by the filled circles, while the ar- 
rows show the derived errors. The best linear fit is indicated as 
the dotted line, while the horizontal soUd line shows the mean 
value. The displayed example is a synthetic simulation. 



B.4. Macroturbulence VmacTo 

The macroturbulenc e con s idered is characterized by a radial- 
tangential model ( iGravl Il975h . In this model a time- 



6 - 



independent velocity field is assumed to be active and modeled 
by assuming that a fraction (A,) of the material is moving radi- 
ally while the complementary fraction moves tangentially (A,). 
In our model for Vmacro, we use the standard assumption that 
the radial and the tangential components are equal in fraction. 
Assuming additionally that the distribution of the velocities in 
the radial and the tangential directions are Gaussian, the effect 
of macroturbulence can be introduced in the synthetic line pro- 
files by means of a convolution with a function, which consists 
of a Gaussian and an additi onal term pr oportional to the error 
function (see equation (4) in lGravlll975h . It changes the shape 
of the line profile, but leaves the equivalent width untouched. 

Only weak lines should be used to determine Vmacro- For 
B-type stars. Si lines are the best suited, since they most ob- 
viously show the presence of Vmacio in their wings. The user 
decides which Si lines should be used throughout the full pro- 
cess. To determine the strength of v^^cm, we convolve the Si 
profiles of the (at that instant) best fitting model with differ- 
ent values of Vmacro^- The consecutive values considered for the 
convolution are chosen using a bisection method, with initial 
steps of 10 kms"' . The bisection continues as long as the step- 
size is larger than or equal to 0.1 kms which will be the 
final precision of the Vmacro determination. For each considered 
I'macro-value, wc computc its log-likelihood in order to quantify 



1 1 .0 
T.f, (kK) 



Figure B.IO. Once a set of possible effective temperatures is 
fixed from the joint predictive power of He I and Si II (see text), 
we can determine for each of these temperatures (in this case: 
10500 and 11 000 K) the correct microturbulence, under the 
assumption that we know the Si abundance. Resulting combi- 
nations are shown for log n(Si)/n(H) = -4.19 (open circles), 
-4.49 (grey filled circles) and -4.79 (black filled circles). The 
displayed example is a synthetic simulation. 



the difference between the obtained synthetic profile and the 
observed line profile, as follows: 



J[-ln(^)-ln(V27?)-i(^)' 



(B.l) 



' At this point the profile has already been convolved with the ap- 
propriate rotational and instrumental profiles. 



with n the considered number of wavelength points within the 
line profile, y, the flux at wavelength point / of the observed 
line profile, yu,- the flux at point / of the synthetic profile (con- 
volved with the Vmacro under consideration), and cr the noise, 
i.e. 1/SNR, of the considered line profile. A thorough discus- 
sion of why this particular likelihood f unction is suitable for 
spectral line fitting can be found, e.g., in lDecin et al.l (l2007h . 
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Figure B.ll. Determination of the macroturbulent velocity of 
P CMa for three different Si lines. The panels should be read 
from left to right as follows. 

Left: The left and right edge of the line region (i.e. outermost 
dots at continuum level) are indicated by the user at the start of 
the procedure. The blue and red edges (grey vertical lines) are 
determined as the point where the flux is half of the minimal in- 
tensity. They determine how far the wings extend towards the 
line center. The black dots represent the points that are finally 
used to determine the macroturbulence. 

Middle: Obtained log-likelihood as a function of Vmacro (see 
Eq. ( IB. 11 1). Note the broad maximum in the log-likelihood dis- 
tribution, which gives rather large error bars. 
Right: The fit with the best Vmacro is shown as the grey profile. 
Note that the significant amount of macroturbulence might be 
expla ined by the fact that th is star exhibits non-radial oscilla- 
tions (iMazumdar et al ]l2006h . 



The same procedure is repeated for all considered line pro- 
files, and the mean, derived from the different lines, gives the 
final Vmacio- The error is set by the maximal deviation of the 
derived Vmacro-values from the final value. From Fig. IB. Ill it is 
clear that the shape of the log-likelihood distribution near the 
maximum can be quite broad, which indicates that the errors in 
Vmacro with rcspcct to the log-likelihood are significant. 

A warning has to be made here. It was recently shown 
by lAerts et al. I (l2009l) that time-dependent line profile varia- 
tions due to non-radial gravity-mode oscillations, which are 
expected in B stars, offer a natural physical explanation for 
the occurrence of macroturbulence. Such oscillations disturb an 
appropriate estimate of v sin / , both from the Fourier method 
and from a goodness-of-fit approach. The only way to assess 
the effect of oscillations on the v sin / determination is by tak- 
ing an appropriate tim e series of line profile variations (e.g., 
Aerts & De Call 2003 ). Luckily, an inappropriate vsin/ esti- 
mate, accompanying a line profile fit in terms of vmacio, does not 
have a serious effect on the determination of the other param- 



eters of the star from snapshot spectra, such that the approach 
implemented in AnalyseBstar is valid as long as the values of 
vsin; and Vmacro are not used for physical interpretations. A 
large value for the macroturbulence is an indication that time- 
dependent pulsational broadening may occur in the star under 
investigation. 



B.5. Surface gravity \ogg 

The surface gravity log^ is the result of fitting the Balmer 
lines. We mainly use Hy and H5, possibly complemented with 
HyS in case of weak winds (i.e. if log Q < -13.80). Since He 
is somewhat blended, and the merging is not always reliable 
for this region (see above), this line will not be used as a pri- 
mary gravity determinator, but will only serve as an additional 
check (together with H/J, in case of stronger winds). Since we 
want to define the profile only from the extreme wings down 
to the strongest curvature (see discussion of the different con- 
tributions below), we excluded the central part of the line, ac- 
counting for several mechanisms, which affect the core of the 
Balmer lines to some extent, e.g., rotation, micro- and macro- 
turbulence and thermal broadening. Once we removed the inner 
part of the line profile, we compare the observed Balmer lines 
with the synthetic line profiles, by considering all possibilities 
for log g which are available at this given effective temperature 
gridpoint. We decide which gravity is the best, by maximiz- 
ing its log-likelihood (see Eq. lB.lb . calculated over all Balmer 
lines. With this procedure, we encountered the following prob- 
lem: Hy and (our main gravity indicators) can sometimes be 
seriously affected by a large amount of blends. Ignoring them 
can lead to a serious overestimation of the gravity, since the 
log-likelihood function will take all these blends into account 
(see Fig. lB.12l i. Therefore, we developed a procedure to find 
the "line continuum", cutting away the blends. The procedure 
is basically a sigma-clipping algorithm which keeps only those 
points that have a flux higher than both neighboring points (the 
so-called 'high points') and where, at the same time, the differ- 
ence is less than 1/SNR. After removing all other flux points, 
only the "line continuum" is left (see middle panel Fig. lB.12l i. 
We interpolate this continuum to obtain the flux at all origi- 
nal wavelength points, and the inner part of the line is added 
again (i.e. the part between the last blend in the blue wing 
and the first blend in the red wing, but obviously still with- 
out the central core). Also all original flux points which devi- 
ate by less than 1/SNR from the interpolated "line continuum" 
are included again. In this way we keep only the flux points 
which really determine the shape of the wing to fit the synthetic 
Balmer wings, and we are able to determine an accurate value 
for log g (see right panel Fig. lB.12l ). We realize that some points 
that may be marked as local continuum in this way, may still 
be lower than what the real local continuum level would be, 
because of the transition of one blend into another However, 
this will only be the case for very few points, which will give 
no significant weight to the log-likelihood. 

The finally accepted surface gravity will be this log^ which 
gives the best fit to all selected Balmer lines simultaneously. 
Half the gridstep in log g could thus be considered a good er- 
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Figure B.12. Illustration of the effect of blends on the deter- 
mination of log g, and the clipping algorithm to improve the 
fit, for p CMa. Left; The many blends in the wings of the pro- 
file prevent an accurate least squares fit and lead to too broad 
wings, implying too high a gravity. Central: From the observed 
line profile (black), only the 'high points' (grey dots), which 
deviate less than 1/SNR from their neighbors, are kept, com- 
plemented with the inner part of the line between the last blend 
in the blue wing and the first blend in the red wing (indicated as 
grey vertical lines), but still excluding the central core, which 
is affected by several broadening mechanisms. Right: After the 
removal of all blends through sigma clipping (see text), we are 
able to obtain a good representation of the line wings (grey 
profile) and, therefore, to derive a very accurate value for the 
gravity. 



ror estimate. However, to account for the coupling with Teff, 
we consider 0. 1 dex as a more appropriate error. This gravity 
should be corrected for centrifugal acceleration (— > \og gc) due 
to stellar rotation by a factor (vsin; )^/Rt, w hen calculating. 



e.g., the mass of a star (iRepolust et al.l 12005 
therein). 



and references 



B.6. Wind parameters 

A change in mass loss rate will mainly affect the shape of Ha. 
For cool objects and weak winds. Ha is nearly 'photospheric' 
and will, in essence, be an absorption profile, with more or less 
symmetric components, since the 're-filling' by the wind is low. 
In case of stronger winds. Ha will take the shape of a typical 
P Cygni profile. For hot objects, where Ha is dominated by 
recombination processes, and high mass-loss rates, the profile 
may even appear in full emission. 

The wind parameter /3 determines the velocity law, which 
directly influences the density. The red wing of Ha is well 
suited to determine /3, since it is formed by emission processes 
alone, averaged over the receding part of the (almost) complete 
wind. For a fixed mass loss, a 'slower' velocity law (i.e., a 
higher /3 value) will result in higher densities in the lower at- 
mosphere, close to the star. This enlarges the number of emit- 
ted photons with velocities close to the line center, resulting in 
more emission. Around the central wavelength, the absorption 
component of the line profile refills and the emission compo- 
nent becomes stronger. Therefore the slope of the red wing of 
the P Cygni profile becomes steeper. In this sense, the steep- 
ness of the red wing is a measure for the value of /3. 
We estimate the wind parameters log Q and /3 by comparing the 
observed Ha profile with the different synthetic Ha profiles 
by making combinations of log Q and /3. We decided to make 
the determination of the wind parameters not too sophisticated, 
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Figure B.13. (Synthetic simulation) For each dif- 
ferent wind strength parameter log Q (7 values, 
from left to right and from top to bottom: log Q - 
-14.30,-14.00, -13.80,-13.60,-13.40,-13.15,-12.70), 
we search for the best /3 (in each panel, the 5 different values 
for /3 are indicated: 0.9, 1.2, 1.5, 2.0, 3.0) by comparing only 
the red wing (grey part of the profile). Then the synthetic 
profiles of each best (log Q, /?) combination are compared to 
the entire Ha profile (bottom right). 



by simply using the best values for log Q and /? available in 
the grid, without interpolation and further refinement. Once we 
have selected for each log Q the best we decide on the best 
combination by considering the best fit to the entire profile. 

6. 7. Some final remarks 

At the end of this iteration cycle, an acceptance test is per- 
formed. When the method can run through the whole cycle 
without needing to update any of the fit parameters, the model 
is accepted as a good model, and gets the flag '2'. If a better 
model was found (i.e., when one or more parameters changed), 
the initial model is rejected and gets flag '0', whereas the im- 
proved model is added to the list of possibilities, and gets flag 
'1' (cf. 'models to check' in the flowchart diagram. Fig. IB. 11 1. 
As long as the list with possibilities with flag ' 1 ' is not empty, 
the parameters of the next possibility are taken as new starting 
values. 

Line profile fits are always performed allowing for a small 
shift in wavelength, or radial velocity (5 wavelength points in 
either direction, corresponding to typically 5kms '), to pre- 
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vent flux differences to add up in case of a small radial velocity 
displacement. 



B.8. The underlying FASTWIND BSTAR06 grid 

Even though the majority of the CoRoT targets consisted of late 
B-type stars, we nevertheless invested in a grid which covers 
the complete parameter space of B-type stars. In this way, we 
hoped to establish a good starting point for future (follow-up) 
detailed spectroscopic analysis of massive stars. The grid has 
been constructed as representative and dense as possible for a 
wide variety of stellar properties within a reasonable compu- 
tation time. Hereafter we give a short summary of the consid- 
ered paramet e rs. For a more detailed description, we refer to 



Lefever et al.1 ( 12007 ah . 



33 effective temperature (Teft): from 10000 K to 32000 K, 
in steps of 500 K below 20000 K and in steps of 1 000 K 
above it. 

on average 28 surface gravities {\ogg) at each effective 
temperature point: from \ogg = 4.5 down to 80% of the 
Eddington limit 

1 'typical' value for the radius (/?») for each (T^ff, \ogg)- 
gridpoint: approximative value from interpolation between 
evolutionary tracks, keeping in mind that a rescaling to the 
'real' value is required when analyzing specific objects. 
Indeed, the real radius can then, in most cases, be deter- 
mined from the visual magnitude, the distance of the star 
and its reddening, and from the theoretical fluxes of the best 
fitting model. 

3 values to fix the chemical composition: n(He)/n(H) = 0. 10, 
0.15, and 0.20, and log n(Si)/n(H) = -4.19, -4.49, and - 
4.79 (i.e. enhanced, solar, and depleted). The background 
elements (responsible, e.g., for radiation pressure and line- 
blanketing) were assumed to have a solar composition. 
7 values for the wind-strength parameter log Q 
1 'typical' value for the terminal wind velocity (vm): esti- 
mated from the relation between ob served terminal velocity 
and photospheric escape velocity (IKudritzki & Pulsl l200d 
and references therein), where the latter quantity has been 
calculated from the actual grid parameters (log^, and 

Teif). 



